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We carry out an extensive numerical study of the dynamics of spiral waves of electrical 
activation, in the presence of periodic deformation (PD) in two-dimensional simulation 
domains, in the biophysically realistic mathematical models of human ventricular tissue 
due to (a) ten-Tusscher and Panfilov (the TP06 model) and (b) ten-Tusscher, Noble, Noble, 
and Panfilov (the TNNP04 model). We first consider simulations in cable-type domains, 
in which we calculate the conduction velocity 6 and the wavelength A. of a plane wave; 
we show that PD leads to a periodic, spatial modulation of 6 and a temporally periodic 
modulation of X; both these modulations depend on the amplitude and frequency of 
the PD. We then examine three types of initial conditions for both TP06 and TNNP04 
models and show that the imposition of PD leads to a rich variety of spatiotemporal 
patterns in the transmembrane potential including states with a single rotating spiral 
(RS) wave, a spiral-turbulence (ST) state with a single meandering spiral, an ST state 
with multiple broken spirals, and a state SA in which all spirals are absorbed at the 
boundaries of our simulation domain. We find, for both TP06 and TNNP04 models, that 
spiral-wave dynamics depends sensitively on the amplitude and frequency of PD and 
the initial condition. We examine how these different types of spiral-wave states can be 
eliminated in the presence of PD by the application of low-amplitude pulses by square- and 
rectangular-mesh suppression techniques. We suggest specific experiments that can test 
the results of our simulations. 
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1. INTRODUCTION 

Sudden cardiac arrest is caused, in many cases, by cardiac arrhyth- 
mias, such as ventricular tachyacardia ( VT) and ventricular fibril- 
lation (VF) (Roger et al., 2011, 2012). Estimates suggest that VF is 
the main reason for death in 30% of the cases in which heart fail- 
ure occurs (Zipes and Wellens, 1998; Fogoros, 2011). Thus, the 
importance of studying such arrhythmias cannot be overempha- 
sized. Such studies must use interdisciplinary approaches because 
they require inputs from biology, bio-medical engineering, and 
cardiology, on the one hand, and physics, non-linear dynam- 
ics, and numerical methods, on the other; methods from these 
areas must be used to study the complicated, non-linear, partial- 
differential-equation models that have been developed for cardiac 
tissue. Such equations can show, inter alia, spiral-wave turbu- 
lence and spatiotemporal chaos, which is believed to be one of 
the mathematical analogs of VF. The study we present here com- 
bines theoretical ideas from spatiotemporal chaos in extended 
dynamical systems with extensive direct numerical simulations, to 
elucidate the effects of periodic deformation (PD) on spiral- wave 
dynamics in detailed mathematical models for cardiac tissue and 
to investigate the elimination of such spiral waves, in the presence 
of PD, by the application of low- amplitude current pulses. 



The mechanisms underlying VT and VF are not understood 
with complete certainty; however, various clinical studies (Fenton 
et al., 2008; Fogoros, 2011) have suggested that such arrhyth- 
mias comprise the abnormal propagation of a wave of electrical 
activation across the ventricles. Such irregular waves appear in 
various ways. For example, they can arise because of an infarc- 
tion scar (De Bakker et al, 1993), which can create an anatomical 
path that anchors electrical waves; they are also seen around 
obstacles in vitro (Valderrabano et al., 2000) and in cell-culture 
experiments (Lim et al., 2006) (anatomically reentry). However, 
reentry can occur in the absence of an anatomical pathway. For 
instance, in vitro experiments (Davidenko et al., 1990; Ikeda et al., 
1996) have shown that fibrillation can be maintained by func- 
tional reentry. However, ectopic activations (Chen et al., 1999; 
Zimmermann and Kalusche, 2001), because of pulmonary veins 
or abnormal cells, can also initiate fibrillation. In particular, both 
experimental (Davidenko et al, 1992; Pertsov et al., 1993; Gray 
et al, 1998; Jalife et al, 1998) and computational (Jalife et al, 
1998; Fenton et al., 2002; Cherry and Fenton, 2008) studies have 
suggested that VT and VF are, respectively, manifestations of (a) 
a rotating spiral (RS) or scroll wave or (b) broken spiral or scroll 
waves that lead to spiral- or scroll-wave turbulence (ST). 
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Several studies have investigated the transition from RS to 
ST, both in experiments on cardiac tissue and in computational 
studies of mathematical models for cardiac tissue; they find that 
this transition can occur because of (a) a steep, increasing ini- 
tial segment in the restitution curve, a plot of the action potential 
duration (APD) versus the diastolic interval (DI) (Roller et al., 
1998; Garfinkel et al, 2000; Fenton et al, 2002), (b) a simi- 
lar steep part in an analogous plot of the conduction velocity 
8 versus DI (Qu et al, 2000b; Fenton et al, 2002), (c) alter- 
nans (Rarma, 1994; Roller et al, 1998; Qu et al., 1999; Cherry 
and Fenton, 2008), and (d) heterogeneities, such as, conduction 
and ionic inhomogeneities (Xie et al., 1998; Shajahan et al., 2007, 
2009; Majumder et al., 2011a,b). Recently, some groups (Zhang 
et al., 2004, 2006; Panfilov et al, 2007; Chen et al, 2008; Weise 
et al., 2011) have begun to study the effects of the deforma- 
tion of cardiac tissue on the RS-ST transition; such a transition 
arises either because of periodic deformation or the stretch- 
activated current associated with such deformation. These studies 
have used simple, two-variable mathematical models for elec- 
trical activation in such tissue. One of the goals of our study 
is to investigate spiral-wave dynamics in general, and RS-ST 
transitions in particular, in a simple mathematical model for 
periodic deformation (PD) of cardiac tissue (Zhang et al, 2004, 
2006; Chen et al., 2008) that we couple with ionically realistic 
human-ventricular-tissue mathematical models, namely, (a) the 
TP06 model, due to ten Tusscher and Panfilov (Ten Tusscher 
and Panfilov, 2006), or (b) the TNNP04 model, of ten Tusscher, 
Noble, Noble, and Panfilov (Ten Tusscher et al., 2004). A defor- 
mation of cardiac tissue leads to modifications of ion-channel 
parameters, because of stretch-activated currents, and also a mod- 
ification of intracellular couplings. Our model for deformation, 
based on Zhang et al. (2004, 2006) and Chen et al. (2008), is 
a very simplified one in which the effects of deformation are 
accounted for only by a temporal modulation of diffusion con- 
stants (see section 2), which are related to intracellular couplings; 
we do not include stretch-activated currents as considered in 
Panfilov et al. (2005, 2007) and Weise et al. (2011). However, 
in spite of this simplified representation of deformation, our 
study yields important results that have been observed in two- 
variable models for cardiac tissue both with periodic deformation 
(PD) (Zhang et al, 2004, 2006; Chen et al, 2008) or mechani- 
cal deformation (Panfilov et al, 2005, 2007; Weise et al, 2011); 
the latter studies include stretch-activated currents. Such stretch- 
activated currents with mechano-electrical feedback have also 
been shown to affect electrical activation in anatomically real- 
istic human models (Reldermann et al., 2010; Ruijpers et al., 
2011); mechano-electrical feedback can enhance electrical acti- 
vation as discussed in Thompson et al. (2011) in the context of 
myofibroblast-myocyte interactions (of course, pharmacological 
and electrochemical interventions can also enhance electrical acti- 
vation). On the positive side, our study uses ionically realistic 
models that have not been employed in such deformation studies 
so far (but see the recent Weise and Panfilov, 2013). We discuss the 
principal results of the studies carried out in Zhang et al. (2004, 
2006), Panfilov et al. (2005, 2007), Chen et al. (2008), and Weise 
et al. (2011), in section 4, where we compare their findings with 
ours. 



We also investigate the efficacy of a low-amplitude suppression 
scheme, developed for the suppression of spiral-wave turbulence 
in 2D models for cardiac tissue (Sinha et al, 2001; Pandit et al, 
2002; Shajahan et al., 2009; Majumder et al, 201 la) in the absence 
of PD. Our study yields several interesting results that we sum- 
marize below before we discuss them in detail. We find first that 
PD leads to a periodic, spatial modulation of the conduction 
velocity 6 and a temporally periodic modulation of the wave- 
length k of a plane wave. We then use three different parameter 
sets, for both TP06 and TNNP04 models, to obtain three dif- 
ferent prototypical spiral configurations, which we use as the 
initial conditions ICl, IC2, and IC3 (see section 2). We find, 
for the TP06 model, that spiral-wave dynamics depends sensi- 
tively on PD and the initial condition. (For similar studies of 
the sensitive dependence of spiral-wave dynamics on inhomo- 
geneities, see Shajahan et al, 2007, 2009.) The initial condition 
ICl can lead to (a) an RS state with n-cycle temporal evolu- 
tion (here n is a positive integer), (b) rotating-spiral states with 
quasi-periodic (QP) temporal evolution, (c) a state with a sin- 
gle meandering spiral MS, which displays spatiotemporal chaos, 
(d) an ST state, with multiple broken spirals, and (e) a quiescent 
state SA, in which all spirals are absorbed; the initial condition 
IC2, with PD, can lead either to (a) an ST state, with multiple 
spirals, or (b) an SA state, with no spirals; and for IC3, it can 
be driven into (a) an ST state, with a single meandering spiral, 
(b) an ST state, with multiple spirals, and (c) the state SA. For 
all these initial conditions, precisely which one of these states is 
obtained depends on the amplitudes A x and A y and the frequen- 
cies f x and f y of the PD in the x and y directions. Spiral-wave 
dynamics in the TNNP04 model, with PD, also shows sensitive 
dependence on PD and the initial condition. This sensitive depen- 
dence on parameters is a hallmark of extended dynamical systems 
that show spatiotemporal chaos (Shajahan et al., 2007, 2009). We 
also study, in the presence of PD, the efficacy of a low-amplitude 
suppression scheme (Sinha et al., 2001; Shajahan et al., 2009) that 
has been suggested, hitherto only without PD, for the suppres- 
sion of spiral-wave turbulence, via low-amplitude current pulses 
applied on a square mesh, in mathematical models for cardiac 
tissue. Furthermore, we develop line-mesh and rectangular-mesh 
variants of this suppression scheme. The latter suppresses spiral 
turbulence in all cases we consider. 

2. METHODS 

The electrical activation of the transmembrane potential V m 
of cardiac tissue is often modeled by a reaction-diffusion-type 
equation, 

dVm lion d 2 Vm d 2 Vm , , 

9f C m dx z 7 dy z 

where C m is the membrane capacitance density, I wn is the sum of 
all the ionic currents that cross the cell membrane, and D x and D y 
are, respectively, the diffusion coefficients along x and y directions; 
such diffusion terms are related to gap junctions (Ten Tusscher 
et al, 2004; Ten Tusscher and Panfilov, 2006), which are net- 
works of protein channels that allow the passage of ions from 
cell to cell. We use two biophysically realistic ionic models for 
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human cardiac myocytes: (a) the ten Tusscher and Panfilov model 
(the TP06 model) (Ten Tusscher and Panfilov, 2006) and (b) the 
ten Tusscher, Noble, Noble, and Panfilov model (the TNNP04 
model) (Ten Tusscher et al., 2004). It is important to check the 
compatibility of the models that we consider in our study with 
the standard version of the original models (Ten Tusscher et al., 
2004; Ten Tusscher and Panfilov, 2006). One way to check for 
such compatibility is by comparing the AP and its morphologi- 
cal properties (Figure SI and Table SI, Supplementary Material 
SI) with those in our studies; we have checked this explicitly, as 
we discuss in detail in the Supplementary Material SI. 

We follow the method suggested in Zhang et al. (2004, 2006) 
and Chen et al. (2008) for the introduction of PD into a mathe- 
matical model for cardiac tissue. In particular, we note that any 
point x = (x, y) in the medium changes to x 7 (f) = (x! (f), / (t)) 
with 



x'(t) =x[l+A x (t)], 
/(t)=y[l+A y (t)], 



(2) 



if we impose a PD with A x (t) = A x cos (2irf x t ) and A y {t) = 
A y cos {Ijtfyt). By substituting Equation (2) into Equation (1), we 
obtain 



, lion 



1 



d z V„ 



[1+A x (t)] 2 X 9x 2 



(3) 



+- 



[1+A/0] 



2^ dy2 ' 



a comparison of Equations (1, 3) shows that Equation (3) can be 
rewritten as 



dV m 

dt 



+ 



D x (t)- 



d z V m 



dx 2 



+ Dy(t)- 



d z V m 



dy 2 



(4) 



withD x (f) = D x (l + A x (t)y 2 andD y (f) = D y (l +A y (t)) 2 . 

In our numerical simulations, we use 2D square domains 
with 1024 x 1024 grid points and lattice spacings Sx = Sy = 
0.25 mm for both TP06 and TNNP04 models, so the sides of 
our square simulation domains are L = 256 mm in the absence 
of PD. We use a forward-Euler method for time evolution, with 
a time step St = 0.02 ms, a five-point stencil for the Laplacian, 
and no-flux (Neumann) boundary conditions. We set the diffu- 
sion coefficients D x = D y = D = 0.00154 cm 2 /ms (Ten Tusscher 
et al, 2004; Ten Tusscher and Panfilov, 2006) for both the TP06 
and the TNNP04 models for our numerical investigations. Other 
parameters for our calculations are given in Tables 1, 2 and an 
examination of the numerical stability of our numerical scheme 
is given in the Supplementary Material SI. 

To examine the spatiotemporal evolution of the electrical sig- 
nal of the transmembrane potential, we obtain the local time 
series of V m (x, y, t), from a representative point (x = 125 mm, 
y = 125 mm) (shown by an asterisk in all pseudocolor plots of 
V m ). To obtain the plots of the inter-beat interval (IBI), we use 
this local time series with 4 x 10 5 data points; the IBI is the time 



Table 1 | Parameters for the periodic deformation (PD) that we use to 
study the wave dynamics in our cable-type and square simulation 
domains in both TP06 and TNNP04 ventricular models. 
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7 



interval between two successive beats in the time series of the sig- 
nal V m . For the power spectra E(co) we use the local time series 
with 2 x 10 5 data points after the initial 10 5 data points have been 
removed to eliminate transients. We present the spatiotemporal 
evolution of V m by a series of (Videos S1-S9) of its pseudocolor 
plots; all these videos use 10 frames per second and each frame is 
separated from the succeeding frame by 8 ms. 

We often have to track the trajectory of the tip of a spiral 
wave in a 2D simulation domain. The tip of such a spiral wave 
is normally defined as the point where the excitation wave front 
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Table 2 | Parameters for the initial spiral-wave configurations for TP06 and TNNP04 models; we refer to these as /CI, /C2, and ICS initial 
conditions. 

Model Initial G Na G kr G ks G pCa G pK a, 

condition (nS/pF) [nS/pF) (nS/pF) (nS/pF) (nS/pF) 

TP06 /C1 14.838 0.153 0.392 0.1238 0.0146 1 

TP06 IC2 5 x 14.838 0.153 0.392 0.1238 0.0146 1 

TP06 IC3 14.838 0.172 0.441 0.8666 0.00219 2 



TNNP04 IC1 14.838 0.096 0.245 0.825 0.0146 1 

TNNP04 IC2 5 x 14.838 0.096 0.245 0.825 0.0146 1 

TNNP04 IC3 3 x 14.838 0.096 0.245 5 x 0.825 0.0146 2 

Here, a f is the scale factor of the time constant tf (Ten Tusscher et al., 2004; Ten Tusscher and Panfilov, 2006). 



t = 600 ms 




* 0 256 512 768 1024 

x (nun) 



FIGURE 1 | Pseudocolor plots of the transmembrane potential V m for 
the TP06 model illustrating plane-wave propagation in a cable-type 
domain, with PD along the axial-direction of the cable, and the 
parameter sets given in Table 1. The Video S1 comprises 21 animations 
that show the spatiotemporal evolution of these plane waves. 



and repolarization wave back meet; this point can be found by a 
variety of methods (Barkley et al., 1990; Fenton and Karma, 1998; 
Fenton et al., 2002; Otani, 2002; Gray et al, 2009; Nayak et al, 
2013). We use the tip-tracking algorithm of Nayak et al. (2013) 
that locates the tip position by monitoring Jnh> the sodium cur- 
rent. Pseudocolor plots of show a fine line along the arm of 
a spiral wave (Figure 2A in Shajahan et al., 2009); this line ter- 
minates in the spiral tip and can, therefore, be used to obtain the 
spatiotemporal evolution of this tip. 

In Figures S2A-F in the Supplementary Material SI, we show 
schematic diagrams for illustrative periodic deformations of a 
small part of our simulation domain, with 5x5 grid points 
(indicated by gpts on the axes of figures); in these diagrams, 
blue, open circles and blue, dashed lines show, at a particular 
instant of time, the deformed simulation domain superimposed 
on the undeformed one, which is represented by black, solid cir- 
cles and black, full lines. We give representative diagrams for the 
case of expansion, with deformations along only x, only y, or 
both x and y directions, in Figures S2A-C, in the Supplementary 
Material SI, at time t = 20ms; the corresponding plots for con- 
traction, at time t = 180 ms, are shown in Figures S2D-F, in the 
Supplementary Material SI. 

3. RESULTS 

We begin by exploring the effects of PD both on plane-wave prop- 
agation and on spiral-wave dynamics; here we vary the oscillation 
amplitude and the frequency in the ranges 0 < A x , A y < 0.5 and 
0 Hz < / < 7.0 Hz; the deformation amplitudes we use are com- 
parable to those in other computational (Zhang et al, 2004, 2006; 
Weise et al, 2011) and experimental (McCULLOCH et al, 1987; 
Noble, 2002) studies; to set the scale of frequencies, we note that 
the frequency of rotation of a single spiral wave is 4.75 Hz for the 
TP06 model and 3.75 Hz for the TNNP04 model (see section 3.2). 
We then study the effects of PD on the suppression scheme of 
Pandit et al. (2002) and Shajahan et al. (2009). 

3.1. PLANE-WAVE DYNAMICS IN A CABLE 

We study plane-wave propagation in a thin, cable-type simulation 
domain, with 16 x 4096 grid points, i.e., L x = 4 mm and L y = 
1024 mm. We inject a stimulus of strength Istimulus = 150pA/pF 
at the left end of the cable for 3 ms and then study the effects of 
PD on the plane wave that propagates through this cable; in par- 
ticular, we measure the conduction velocity 6 and wavelength X 



of the propagating wave in the cable. We find that 0 ~ 70.6 cm/s 
and X ~ 21.6 cm for the TP06 model, and 8 ~ 67.8 cm/s and 
X ~ 18.9 cm for the TNNP04 model, in the absence of PD. As sug- 
gested in Clayton and Panfilov (2008), Shajahan et al. (2009), and 
Ten Tusscher et al. (2004), it is useful to test the accuracy of the 
numerical scheme by varying both the time and space steps that 
we use for integration. We illustrate this for the TP06 model by 
measuring 9 for a plane wave, which is injected into the medium 
by stimulating the left boundary of our simulation domain. We 
find that, with Sx = 0.025 cm, 6 increases by 1.6% as we decrease 
St from 0.02 to 0.01 ms; if we use St = 0.02 ms and decrease Sx 
from 0.025 to 0.015 cm then 6 increases by 4.7%; such changes 
are comparable to those found in earlier studies (Ten Tusscher 
et al, 2004; Shajahan et al, 2009). 

In Figure I(a00-a20) we show, at time t = 600 ms, when 
PD is applied along the axial direction of the cable, pseudo- 
color plots of the transmembrane potential V m for the TP06 
model with PD along the axial direction of the cable, and the 
parameter sets given in Table 1. The Video SI comprises 21 
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animations that show the spatiotemporal evolution of the plane 
waves in Figure I(a00-a20); these animations and Figure l(a00- 
a20) show that the conduction velocity 9 is modulated in space 
and the wavelength k is modulated in time because of the PD. 
Figure 2 illustrates these modulations via plots of 9p and 9b ver- 
sus x for the conduction velocities of the wave front (Figure 2A) 
and the wave back (Figure 2B), respectively; here the subscripts 
F and B stand for wave front and wave back, respectively; and 
Figure 2C shows the corresponding plot for k versus time t; 
in these plots we use the representative PD parameter values 
A x = 0.3 and f x = 5.0 Hz for the TP06 model. We calculate the 
conduction velocities 9f(x) and 6b(x), in the cable-type domain 
with PD, by recording the positions of the wave front and the 
wave back at times f and t + At, with At = 2 ms; the wave-front 
and wave-back conduction velocities, at the point x at time f, 
are 9f (x) = Afx/ At and 9b(x) = Agx/At, where Afx and Agx 
are, respectively, the distances traveled by the wave front and 
wave back in the time interval Ar. We locate the position of the 
wave front by finding the value of x at which V m ~ 0 mV; we 
define the position of the wave back as the point, behind the 
wave front, at which a secondary action potential can just be 
initiated by an additional stimulus (this turns out to occur at 



a value of V m that is ~ 75% of the repolarization phase of the 
action potential). We obtain the wavelength k(t) by measuring 
the distance between the wave front and the wave back at 
time t. 

In Figures 2A,B, the open circles show the values of 9p(x) and 
9b (x) , respectively, that we obtain by the method described above; 
the red lines show smooth sinusoidal envelopes, with amplitude 
~ 31.2 cm/s and spatial period ~ 14.5 cm, that give the average 
modulations of these conduction velocities with x. Note that, in 
the absence of PD, 9 ~ 70.6 cm/s (this is shown via a gray, dashed 
line in Figures 2A,B); therefore, the electrical wave can travel 
~ 70.6// cm in l//s; hence, for a given PD frequency/, the spa- 
tial period of oscillation of 9f(x) and 9b(x) is ~ 70.6// cm; the 
representative plots of Figures 2A,B, in which / = 5 Hz and the 
period is ~ 70.6/5 = 14.12 cm, are consistent with this estimate. 

Figure 2C shows that k is a periodic function of t with a period 
t; we expect that r = 1//, where / is the PD frequency; the illus- 
trative plot in Figure 2C, with / = 5 Hz, is consistent with this 
expectation because r ~ 202 ms; the gray, dashed line shows the 
value of k that we obtain in the absence of PD. 

It is useful to study how 9 and k of a plane wave behave, 
in the presence of PD, when we change the values of the time 
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FIGURE 2 | The spatial modulation of 9 and the temporal modulation of 
X of a plane wave propagating in a cable-type domain with PD. Plots 
versus distance x along the cable of the conduction velocities with 
D= 0.001 54 cm 2 /ms, St = 0.02 ms, and Sx = 0.25 mm (A) 0 F , of the wave 
front, and (B) 0b, of the wave back; (C) plot versus time t of the wavelength 
X. The exact analogs of (A-C) are shown in (D-F), for D = 0.001 54 cm 2 /ms, 
St = 0.01 ms, and Sx = 0.25 mm, and in (G-l), for D = 0.00077 cm 2 /ms, 



St = 0.02 ms, and Sx = 0.25 mm. We use the representative PD parameter 
values A x = 0.3 and t x = 5.0 Hz for the TP06 model. Open circles show the 
values from our calculation; the red lines show smooth sinusoidal envelopes; 
in the absence of PD, 0 ~ 70.6 cm/s, ~ 71 .7 cm/s, and ~ 47 cm/s, 
respectively, for above three parameter sets [gray, dashed lines in (A,B)I; in 
(C) the gray, dashed line shows the value of X that we obtain in the absence 
of PD. 
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step and the diffusion coefficients. We find that, in the pres- 
ence of PD, 9 and k continue to oscillate, as in Figures 2D-I, 
as functions of x and f, respectively; the mean values of 9 and 
X, about which these oscillations occur, are close to their val- 
ues without PD, which depend on the diffusion coefficients and 
marginally on the time step: In Figures 2D,E we show, for D = 
0.00154 cm 2 /ms, St = 0.01 ms, and Sx = 0.25 mm, the analogs of 
Figures 2A-C; and in Figures 2G-I we give their counterparts for 
D = 0.00077 cm 2 /ms, St = 0.01 ms, and Sx = 0.25 mm. 

The TNNP04-model analogs of the TP06-model 
Figure I(a00-a20) are given in Figure S4(a00-all) in the 
Supplementary Material SI. Our results for the TNNP04 model 
are similar to those for the TP06 model. 

3.2. SPIRAL-WAVE DYNAMICS IN A HOMOGENEOUS DOMAIN 

We move now to systematic studies of spiral-wave dynamics in 
a 2D, square simulation domain with side L = 256 mm, in the 
presence of PD, for both TP06 and TNNP04 models. 

In the absence of PD, two methods are used to initiate spi- 
ral waves in simulations (Pertsov et al., 1993; Bernus et al., 2002; 
Ten Tusscher et al., 2004; Shajahan et al, 2009) and experiments 
(Davidenko et al, 1992; Pertsov et al, 1993), namely, (1) the SI, 
S2 cross-field protocol and (2) the SI, S2 parallel-field protocol. 



We describe in the Supplementary Material SI the precise SI, 
S2 cross-field protocol that we use to obtain spiral waves in our 
simulation domain. We use three types of spiral-wave initial con- 
figurations for our subsequent studies; we refer to these as ICl, 
ICl, and 7C3 initial conditions (see Table 2 for parameter val- 
ues). In Figures S5A-C, in the Supplementary Material SI, we 
show the time evolution of pseudocolor plots of V m for the TP06 
model with the ICl initial configuration; similar plots are shown 
in Figures S5D-E and (F,I), respectively, for the ICl and IC3 ini- 
tial configurations. The TNNP04-model analogs of Figures S5A-I 
are given in Figures S6A-I, in the Supplementary Material SI. 

In Figures 3A-C, we show pseudocolor plots of V m at times 
t = 0 s, f = 2 s, and t = 4 s, respectively, for the initial condition 
ICl in the TP06 model, in the absence of PD; this initial con- 
figuration evolves to a state with a rotating spiral (RS) in the 
medium; the animation (a) in Video S2 shows the spatiotemporal 
evolution of V m for this case. The local time series of V m (x, y, t), 
from the representative point (x = 125 mm, y = 125 mm) (the 
asterisk in Figure 3C), is shown in Figure 3D for 2 s < t < 6 s; 
a plot of the IBI is given in Figure 3E, which shows that, after 
initial transients (roughly the first 10 beats), the spiral wave 
rotates periodically with an average rotation period T — 210 ms. 
In Figure 3F, we plot the power spectrum E(w), which we have 






FIGURE 3 | Spatiotemporal evolution of V m for the initial condition /CI 
for the TP06 model in the absence of PD. (A-C) Pseudocolor plots of V m at 
times f = 0 s, f = 2 s, and f = 4s, respectively, showing the evolution toward 
a state with a rotating spiral (RS); the animation (A) in Video S2 shows the 
spatiotemporal evolution of V m for this case. (D) The local time series of 
V m (x, y, t), from the representative point (x = 125 mm, y = 125 mm) (the 
asterisk in (C)) for 2 s < f < 6s; (E) a plot of the IBI, which we obtain from 



this time series, of length 4 x 10 5 iterations; (F) the power spectrum EM, 
obtained from the local time series of (D), with discrete peaks at the 
fundamental frequency cof ~ 4.75 Hz and its harmonics. The spiral-tip 
trajectory traces a roughly circular path, with radius l c ~ 20 mm, which is 
shown, for 3.6 s < f < 4 s, by the white line that has been superimposed on 
the pseudocolor plot of V m in (C); a magnified view of this path is 
shown in (G). 
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obtained from the local time series of V m mentioned above; dis- 
crete peaks in E(co) appear at the fundamental frequency cof ~ 
4.75 Hz and its harmonics. The periodic nature of the local time 
series of V m , the flattening of the IBI, and the discrete peaks 
in E(w) show that the temporal evolution of the spiral wave is 
periodic; therefore, the spiral-tip trajectory traces a roughly cir- 
cular path with radius l c ~ 20 mm; this circular path is shown, 
for 3.6 s < t < 4 s, by the white line that has been superimposed 
on the pseudocolor plot of V m in Figure 3C; an expanded ver- 
sion of this path is shown in Figure 3G. The Figures S7, S8 in the 
Supplementary Material SI show the analogs of Figure 3 for ini- 
tial conditions IC2 and Id, respectively; and the animations (b) 
and (c) in Video S2 show the spatiotemporal evolution of V m for 
these cases. These animations, the pseudocolor plots of V m , the 
representative local time series of V m , the plots of the IBI, and the 
power spectra show that the initial conditions IC2 and IC3 lead, 
respectively, to spatiotemporal chaos and spiral turbulence (ST), 
with a single spiral meandering chaotically, and broken spirals, 
respectively, in the simulation domain. 

Figures S9A-G, S10A-G, and S11A-F (Supplementary 
Material SI) show, respectively, the TNNP04 analogs of the TP06 
Figures 3A-G (for /CI), Figures S7A-G in the Supplementary 
Material SI (for IC2), and Figures S8A-F in the Supplementary 
Material SI (for 7C3); the spatiotemporal evolution of V m for 
these three initial conditions for the TNNP04 model are given in 
animations (d), (e), and (f) in Video S2. From these animations 
and the Figures S9A-G, S10A-G, S11A-F (Supplementary 
Material SI) we conclude that the spatiotemporal evolution 
of V m in the TNNP04 model, without PD, is similar to, but 
not identically the same as, that in the TP06 model for the 
initial conditions ICl, IC2, and JC3. One difference is that, in 
the TNNP04 model, we have a Z-type, spiral-tip trajectory in 
Figures S10C,G (Supplementary Material SI), whereas, for the 
same initial condition, we have an open spiral-tip trajectory 
(Figures S7C,G in the Supplementary Material SI) in the TP06 
model. This shows that spiral-wave dynamics in these two 
models, without PD, depends sensitively on the ionic details of 
these models and on the initial conditions. 

3.3. SPIRAL WAVES WITH PD 

We present systematic studies of spiral-wave dynamics here by 
using ICl, IC2, and IC3 initial configurations in the presence 
of PD. In the TP06 model, these initial configurations lead, 
respectively, to (a) an RS state with a roughly circular spiral-tip 
trajectory, (b) a single meandering spiral with turbulence (we 
refer to this as SMST henceforth), and (c) multiple-spiral tur- 
bulence (MST) with broken spiral waves in the absence of PD, 
as we have described above. For the TNNP04 model the analogs 
of these states are (a) RSC, a state with a rotating spiral whose 
tip trajectory is roughly circular, (b) RSZ, a state with a rotat- 
ing spiral whose tip trajectory is roughly Z-type, and (c) an MST 
state. 

We first consider the time evolution of ICl for the TP06 
model in the presence of PD, for which we deform the medium 
periodically along both x and y directions, with amplitudes and 
frequencies in the ranges 0.1 < A x = A y < 0.5 and 1.0 Hz < 
fx =f y < 7.0 Hz, respectively. 



In Figures 4A-D we show pseudocolor plots of V m at time 
t = 4s for (a) A x = A y = 0.1, f x =f y = 1.0 Hz, (b) A x = A y = 
0.1, f x =f y = 3.0Hz, (c) A x = A y = 0.1, f x =fy = 5.0Hz, and 
(d) A x = Ay = 0.1, f x =f y = 7.0 Hz, respectively. The RS state, 
which we obtain in the absence of PD, does not evolve into an 
MST state in cases (a), (b) and (c); however, in case (d) the spi- 
ral arm splits into multiple spirals to yield an MST state with 
mild spatiotemporal chaos, in so far as the dominant spiral does 
not break down but continues to evolve somewhat like a mother 
rotor (Sarnie and Jalife, 2001; Chen et al., 2003; Wu et al., 2004; 
Ideker and Rogers, 2006); the animations (al), (bl), (cl), and 
(dl) in Video S3 show the spatiotemporal evolution of V m for 
these cases in the time interval 0s < f < 4s; this video uses 10 
frames per second (fps) and each pseudocolor plot of V m is sep- 
arated from its predecessor by 8 ms. The spiral-tip trajectories, 
which follow from this spatiotemporal evolution, are shown in 
Figures 3E-H for 3.6 s < t < 4 s (in Figure 4H we give the tip 
trajectory for the main, central spiral in Figures 4A-D); these 
tip trajectories are nearly circular with radii l c ~ 18 mm, but, 
as we show below, the temporal evolution of V m is different in 
these cases. To examine this evolution, we obtain the local time 
series of V m (x, y, t), from the representative point (x = 125 mm, 
y = 125 mm) (the asterisks in Figures 4A-D), and therefrom the 
plots of the IBI shown in Figures 4I-L and the power spectra 
of Figures 4M-P. Discrete peaks in E(a>) appear at the funda- 
mental frequency u>f ~ 4.75 Hz and a few other frequencies (see 
the caption of Figure 4). From Figures 4I-L we see that the IBI 
displays a slight upward trend; this implies that, although the 
temporal evolution is nearly periodic, there is a gentle drift, 
toward lower frequencies, in the rotation rate of the dominant 
spiral. Furthermore, there are small oscillations in the IBI in 
Figure 41 (a 5-cycle), (j)(a 3-cycle), and (l)(a 2-cycle), but not in 
Figure 4K (a 1-cycle); the natures of these oscillations and their 
cycle lengths are confirmed by the Poincare-type return maps, 
shown in Figures 4Q-T and corresponding to the IBI plots in 
Figures 4I-L, respectively; in these return maps, successive points 
are connected by lines. 

Similarly, we study the dependence of spiral-wave dynamics on 
the amplitudes A x and A y of the PD, with the frequencies f x =f y 
held at a fixed value. In Figure 5 we show the pseudocolor plots 
ofV m attimef = 4sfor(a)A x = A y = 0.1, f x = f y = 1.0Hz,(b) 
A x =A y = 0.2, fx = f y = 1.0 Hz, (c) A x =A y = 0.3, f x =f y = 
1 .0 Hz, and (d) A x = A y = 0.4, f x = f y = 1.0 Hz. The spiral wave 
does not split into multiple spirals for these representative val- 
ues of the amplitudes and frequencies. The animations (al), (a2), 
(a3), and (a4) in Video S3 show the spatiotemporal evolution of 
these spiral waves for the interval 0 s < t < 4 s. To examine this 
evolution, we obtain the local time series of V m {x, y, t), from the 
representative point (x = 125 mm, y = 125 mm) (the asterisks 
in Figures 5A-D) and the corresponding tip trajectories of spiral 
waves, in the time interval 3.6s<f<4s (blue lines with black 
points in Figures 5E,H, respectively); these tip trajectories trace 
nearly circular paths, with radii l c — 18 mm in Figures 5E,F; they 
are of the meandering type in Figures 5G,H, with linear extents 
l c ~ 24 mm and l c ~ 75 mm, respectively. From the local time 
series of V m mentioned above, we obtain the plots of the IBI 
shown in Figures 4I-L and the power spectra of Figures 4M-P; 
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FIGURE 4 | Time evolution of the RS state in the TP06 model in the 
presence of PD with a fixed amplitude. Pseudocolor plots of V m at time 
f = 4s for (A) A x = A y = 0.1, f x = f Y = 1.0 Hz, (B) A x = A Y = 0.1, f x = 
f K = 3.0Hz, (C) A x = A Y = 0A, f x = f y = 5.0 Hz, and (D) A x = A y = 0A, 
f x = f y = 7.0 Hz, respectively; the animations (a1), (b1), (c1), and (d1) in 
Video S3 show the spatiotemporal evolution of V m for these cases in the 
time interval 0s< f <4s. (E-H) Spiral-tip trajectories, which follow from 
these spatiotemporal evolutions, for 3.6s < f < 4s [in (H) we give the tip 
trajectory for the main, central spiral (D)]. We obtain the local time series 
of V m (x, y, t), from the representative point (x= 125 mm, y = 125 mm) 
[the asterisks in (A-D)[, and therefrom the plots of the IBI (l-L) and the 
power spectra (M-P). The discrete peaks in EM appear at the following 



frequencies: (M) &>-, = 4.75 Hz, a> 2 = 9.5 Hz, to 2 = 14.25 Hz, (N) 
to! = 4.75 Hz, (o 2 = 9.5 Hz, m 2 = 14.25 Hz, and small peaks at a>i = 3 Hz, 
co2 = 7.75 Hz, m = 11 Hz, a 4 = 12.5 Hz, w b = 15.75 Hz, (O) mi = 4.75 Hz, 
m 2 = 9.5 Hz, a> 2 = 14.25 Hz, and (P) coy = 4.75 Hz, to 2 = 9.5 Hz, 
cu2 = 14.25 Hz. In (l-L) we see that the IBI shows a slight upward trend; 
this implies that, although the temporal evolution is nearly periodic, there 
is a slight drift, toward lower frequencies, in the rotation rate of the 
dominant spiral; note also the mild oscillations in the IBI in (I) a 5-cycle, 
(J) a 3-cycle, and (L) a 2-cycle, but not in (K) a 1 -cycle; the natures of 
these oscillations and their cycle lengths are confirmed by the 
Poincare-type return maps, shown in (Q-T), respectively; in these return 
maps, successive points are connected by lines. 



discrete peaks in E(a>) appear at the fundamental frequency cot — 
4.75 Hz and at the frequencies listed in the caption of Figure 5; 
these peaks indicate that, in Figures 5M,N, we also have some 
high-order cycles; the broad-band power spectra in Figures 50,P 



provide evidence for spiral turbulence with a meandering spiral 
(SMST). In Figures 5Q-T, we show Poincare-type return maps 
which we obtain from the IBI plots in Figures 5I-L; in these 
maps successive points are connected by lines. These plots give 
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FIGURE 5 | Time evolution of the RS state in the TP06 model in the 
presence of PD with a fixed frequency. Pseudocolor plots of V m at 
time f = 4 s for (A) A x = A y = 0.1, f x = f y = 1.0 Hz, (B) A x = A Y = 0.2, 
f x = f y = 1 .0 Hz, (C) A x = A y = 0.3, f x = f y = 1 .0 Hz, and (D) A x = A Y = 
0.4 , f x = fy = 1.0 Hz. The animations (a1-a4) in Video S3 show the 
spatiotemporal evolution of these spiral waves for the interval 
0s<f<4s. (E-H) Spiral-tip trajectories, which follow from these 
spatiotemporal evolutions, for 3.6s<f<4s. We obtain the local time 
series of V m (x, y, f), from the representative point (x= 125 mm, 



y= 125mm) (the asterisks in Figures 4A-D) and therefrom the plots of 
the IBI shown in (l-L), the power spectra EM in (M-P), and the 
Poincare-type return maps, which we obtain from the IBI plots and 
which show (Q) a 5— cycle, (R) a 5— cycle, (S) chaotic behavior, and (T) 
chaotic evolution; discrete peaks in EM appear at the following 
frequencies: (M) a>, =4.75 Hz, a> 2 = 9.5Hz, tu 2 = 14.25 Hz, (N) 
a>-|=4.75Hz, m 2 = 9.5 Hz, to 2 = 14.25 Hz and small peaks at 
(u 1= 3.75Hz, m 2 = 8.5 Hz, to 3 = 13.25 Hz, o) 4 = 15Hz, tu 5 = 17.75Hz, (O) 
a> 1= 4.75Hz, a> 2 = 9.5 Hz, and (P) ^ = 4.5 Hz, to 2 = 9-5Hz. 



additional evidence for five cycles in Figures 5I,J,M,N and of 
chaotic behavior in Figures 5K,L,0,P. The lines in Figures 5Q,R 
move from the bottom-left corner to the top-right corner; this 
suggests a low-frequency modulation of the spiral-wave dynam- 
ics because of the PD; this is associated with the upward trend in 
the IBI plots of Figures 5IJ. 



We focus next on the types of ST states that we obtain, with 
PD applied along both x and y axes, when we start with the 
ICl initial condition. In Figure 6 we show three representative 
ST states; Figures 6A-C show, respectively, pseudocolor plots of 
the transmembrane potential V m for PD with (a) A x = A y = 
0.3, f x =f y = 3.0 Hz, (b) A x = A y = 0.3, f x = f y = 5.0 Hz, and 
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FIGURE 6 | Temporal evolution of representative ST states in the TP06 
model with PD along both spatial directions. Three ST states, which we 
obtain with the initial condition /CI and PD, are shown via pseudocolor plots of 



the transmembrane potential V m with (A) A x 



A, : 



: 0.3, f x =f y = 5.0 Hz, and (C) A x 



-- 0.3, f x 
■ 0.4, f x 



= 3.0 Hz, (B) 

■ 7.0 Hz; the 



animations (b3), (c3), and (d4) in Video S3 show, respectively, the 
spatiotemporal evolution of V m for these cases in the time interval 0 s < f < 4s. 
We obtain the local time series of V m (x, y, f), from the representative points 
(x = 125 mm, y = 125 mm) and (x = 50 mm, y = 50 mm), shown by asterisks 



in (A-C); from these local time series, we obtain the plots of the IBI (D-F) and 
the power spectra (G-l), with open-blue and black-filled circles for the time 
series from (x = 125 mm, y = 125 mm) and (x = 50 mm, ^ = 50 mm), 
respectively. These pseudocolor plots and animations of V m and the plots of the 
IBI and power spectra show that we have, roughly speaking, three types of ST 
states with (A) multiple spirals, (B) a stable spiral core with broken spiral arms, 
and (C) a single dominant meandering spiral; the second case (B) displays a 
coexistence of a quasiperiodic and an ST state because of the dominant spiral 
at the center and the broken spirals generated from its arm. 



(c) A x = Ay = 0.4, f x =fy = 7.0 Hz; and the animations (b3), 
(c3), and (d4) in Video S3 show, respectively, the spatiotempo- 
ral evolution of V m for these cases in the time interval 0 s < 
t < 4 s. To examine this evolution, we obtain the local time 
series of V m (x, y, t ), from the representative points (x = 125 mm, 
y = 125 mm) and (x = 50 mm, y = 50 mm), both of which are 
indicated by asterisks in Figures 6A-C; from these local time 
series, we obtain the plots of the IBI shown in Figures 6D-F and 
the power spectra of Figures 6G-I, with open-blue and black- 
filled circles for the time series from (x = 125 mm, y = 125 mm) 
and (x = 50 mm, y = 50 mm), respectively. These pseudocolor 
plots and animations of V m and the plots of the IBI and power 
spectra show that we have, roughly speaking, three types of 
ST states with (a) multiple spirals (Figure 6A), (b) a stable 
spiral core with broken spiral arms (Figure 6B), and (c) a sin- 
gle, dominant, meandering spiral (Figure 6C); the second case 
(b) displays a coexistence of a quasiperiodic and an ST state 
because of the dominant spiral at the center and the broken 
spirals generated from its arm. Such coexistence behaviors have 
been observed in both computational (Xie et al., 2001; Fenton 
et al., 2002; Cherry and Fenton, 2008) and experimental studies 



(Nash et al., 2006; Masse et al., 2007), which include in vivo 
experiments. 

We also obtain quiescent (Q) states with no spirals because 
of the absorption of spiral waves at the boundaries as shown in 
the pseudocolor plots of V m in Figure S12 in the Supplementary 
Material SI and in Video S3. 

To illustrate the rich variety of spatiotemporal patterns, we 
summarize our results for the TP06 model, with the initial con- 
dition 7C1, by presenting a selection of pseudocolor plots of 
V m in Figure S13(al-d5) in the Supplementary Material SI (for 
parameter sets see Table 1). The animations in Video S3 show the 
spatiotemporal evolution of V m for these cases in the time interval 
0 s < t < 4 s. To examine this evolution, we obtain the local time 
series of V m (x, y, t), from the representative points (x = 125 mm, 
y= 125 mm); these are shown in Figure S14 (Supplementary 
Material SI); from these local time series, we obtain the plots of 
the IBI (Figure S15 in the Supplementary Material SI) and the 
power spectra (Figure S16 in the Supplementary Material SI). 

The counterparts of Figures S13-S16 in the Supplementary 
Material SI, for initial conditions IC2 for 7C3, are given, respec- 
tively, in Figures S17-S24 in the Supplementary Material SI. 
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For the initial conditions IC2 and IC3 the analogs of the ani- 
mations in Video S3 are given, respectively, in Videos S4, S5. For 
IC2, with PD along both axes and different values of the ampli- 
tude and the frequency, we examine the time series of V m (x, y, t), 
from a representative point in the simulation domain (Figure S18 
in the Supplementary Material SI), the plots of the IBI (Figure 
S19 in the Supplementary Material SI) and the power spec- 
trum E(w) (Figure S20 in the Supplementary Material SI), and 
the spatio temporal evolution of V m (given by the animations in 
Video S4) and conclude therefrom that, in this case, we obtain 
either (a) a Q state with no spirals [see animations (a4), (a5), (b4), 
(b5), (d4) and (d5) in Video S4] or (b) an MST state with broken 
spiral waves (see the remaining animations in Video S4). A sim- 
ilar analysis, for IC3 and PD along both x and y axes, based on 
time series of V m (Figure S22 in the Supplementary Material SI), 
plots of the IBI (Figure S23 in the Supplementary Material SI), 
the power spectrum (Figure S24 in the Supplementary Material 
SI), and the spatiotemporal evolution of V m (the animations in 
Video S5) suggests that here we can have (a) a Q state with no 
spirals [see animations (b3) (b4), (b5), (c4), (d4), and (d5) in 
Video S5], (b) an SMST state [see animations (al) and (a5) in 
Video S5], or (c) an MST state with broken spiral waves (see the 
rest of the animations in Video S5). 

The TNNP04 model with PD also exhibits a rich variety of 
spatiotemporal patterns with spiral waves like the TP06 model 
as we discuss in the Supplementary Material SI (see Figures 
S25-S36). These figures and the associated Videos S6-S8 show 
that spiral-wave dynamics with PD in the TNNP04 model is 
quantitatively different from, but qualitatively similar to, that in 
the TP06 model. These differences arise because of the differences 
in the calcium-ion dynamics in these models; such dynamics can 
play an important role in the spatiotemporal evolution of spiral 
waves (Weiss et al., 2005; Ter Keurs and Boyden, 2007). 

We have discussed spiral- wave dynamics in TP06 and TNNP04 
models in the presence of PD along both x and y directions, with 
initial conditions of types ICl, IC2, and IC3. We have also carried 
out systematic simulations of spiral- wave dynamics in both these 
models, with PD along only one (say x) direction. Here too our 
results are, in the main, qualitatively similar to those we have pre- 
sented above. Of course, there is anisotropic diffusion if the PD is 
only along one direction. However, Q, RS, and ST states appear; 
an overview of their spatiotemporal evolution is given in Figures 
S37-S42 in the Supplementary Material SI. 

3.4. SUPPRESSION OF SPIRAL WAVES 

One of the principal goals of our extensive numerical studies 
of spiral-wave dynamics in the TP06 and TNNP04 models with 
PD is to understand its role in enhancing or suppressing spiral- 
wave turbulence; this is an important step in developing an 
effective, low-amplitude suppression technique for the elimina- 
tion of turbulence with single or multiple spirals. So far, various 
low-amplitude suppression algorithms have been developed to 
eliminate spiral waves in monodomain mathematical models of 
cardiac tissue; in these algorithms the control pulses are applied 
in several ways. These include the following: (a) periodic stim- 
ulation at a point (Zhang et al., 2003; Yuan et al, 2005); (b) a 
line stimulus that must be applied to one of the boundaries (Tang 



et al., 2008; Miguel et al., 2009); (c) an array of low-voltage control 
pulses, which must be swept over the simulation domain (Sinha 
and Sridhar, 2007; Sridhar and Sinha, 2008); or (d) the mesh- 
based, low-amplitude suppression scheme we describe (Sinha 
et al, 2001; Pandit et al., 2002). We have provided an overview of 
such low-amplitude suppression schemes, in the absence of PD, in 
earlier studies (Shajahan et al., 2009; Nayak, 2013); the most suc- 
cessful of these is based on a mesh-based suppression algorithm; 
this suppression scheme (Shajahan et al., 2009; Majumder et al., 
2011a; Nayak et al., 2013) can suppress spiral waves of electrical 
activation even in the presence of conduction, ionic, and fibrob- 
last heterogeneities (Shajahan et al., 2009; Majumder et al., 201 la; 
Nayak et al, 2013). We now investigate the efficacy of this mesh- 
based suppression scheme for both TP06 and TNNP04 models in 
the presence of PD. 

In this mesh-based suppression scheme, we apply a current 
pulse of amplitude 75 pA/pF for 0.2 s over a mesh that divides our 
square simulation domain with L = 256 mm into 64 square cells 
of side I = 32 mm each; this pulse makes the links of the mesh 
refractory and, thereby, effectively imposes Neumann boundary 
conditions for any block inside the mesh; therefore, spiral waves 
inside a block are absorbed on the links of the mesh that bound 
the block. We have also extended this mesh-based scheme to one 
that uses control pulses on a set of parallel lines; in this line- 
based scheme, we apply a current pulse of amplitude 125 pA/pF 
for 0.6 s over a set of parallel lines separated from each other 
by I = 32 mm. As we show in the Supplementary Material SI 
(see Figures S43, S44), both these schemes succeed in suppressing 
spiral-wave turbulence in the TP06 and TNNP04 models with- 
out PD; the line-based scheme uses a higher amplitude for the 
control pulse and a longer duration of application than the mesh- 
based one because the former has fewer control-pulse segments 
than the latter. Furthermore, we show (see Figures S45-S48, in 
the Supplementary Material SI) that the line-based scheme works 
with PD only if the PD is applied along one spatial direction. 

A minor modification of our line-based suppression scheme 
suppresses spiral-wave turbulence: we use a rectangular-mesh- 
based control scheme, in which we add a few control lines 
perpendicular to the parallel lines of the line-based suppression 
scheme. We present a comparison of spiral-wave suppression by 
low-amplitude pulses on square, line, and rectangular suppres- 
sion meshes in the TP06 model, with PD along both x and y 
directions: We impose PD along both x and y directions with the 
illustrative amplitudes A x = A y = 0.3 and frequencies/^ = f y = 
5 Hz for the initial configurations ICl, IC2, and IC3 (pseudo- 
color plots of V m in Figures 7A,E,I, respectively). We apply the 
following control pulses: (1) pulses with amplitude 75pA/pF for 
t = 0.2 s over a square mesh (Figures 7B,F,J), with each square 
block of side I = 32 mm; (2) pulses with amplitude 125pA/pF 
for t = 0.6 s over a line mesh (Figures 7C,G,K), with inter-line 
spacing I = 32 mm; (3) pulses with amplitude 125 pA/pF for t = 
0.6 s over a rectangular mesh (Figures 7D,H,L), with block sides 
l x = 32 mm and l y = 64 mm. These pseudocolor plots of V m and 
the associated animations in Video S9 show that such spiral-wave 
states, with ICl, IC2, and IC3 initial conditions, are suppressed by 
both square- and rectangular-mesh suppression but not by line- 
mesh suppression. Our rectangular-mesh suppression scheme is a 
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FIGURE 7 | Comparison of spiral-wave suppression by low-amplitude 
pulses on square, line, and rectangular control meshes in the TP06 
model, with PD along both xand y directions. We impose PD along both 
x and y directions with the illustrative amplitudes A x = A y = 0.3 and 
frequencies f x = f y = 5 Hz for the initial configurations /CI , IC2, and IC3 
(pseudocolor plots of V m in (A,E,I), respectively). We apply the following 
control pulses: amplitude 75 pA/pF for f = 0.2 s over a square mesh (B,F,J), 
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with each square block of side / = 32 mm; amplitude 125 pA/pF for f = 0.6 s 
over a line mesh (C,G,K), with inter-line spacing /= 32 mm; amplitude 
125 pA/pF for f = 0.6 s over a rectangular mesh (D,H,L), with block sides 
l x = 32 mm and l y = 64 mm. These pseudocolor plots of V m and the 
associated animations in Video S9 show that these spiral states, with /CI , 
IC2, and IC3 initial conditions, are suppressed by both square- and 
rectangular-mesh control but not line-mesh control. 



significant improvement over the square-mesh one because it uses 
fewer control lines than the latter. The results of similar studies 
for the TNNP04 model are given in Figure S49 (Supplementary 
Material SI). 

4. DISCUSSION AND CONCLUSION 

We have carried out detailed and systematic numerical studies of 
the effects of periodic deformation (PD) on spiral-wave dynam- 
ics in ionically realistic mathematical models for cardiac tissue 
by introducing PD in the recently developed TP06 and TNNP04 
mathematical models for human ventricular tissue (Ten Tusscher 
et al., 2004; Ten Tusscher and Panfilov, 2006), in which we include 
deformation as in Zhang et al. (2004, 2006) and Chen et al. 
(2008). We also investigate, in 2D simulations with PD, the effica- 
cies of square-, rectangular-, and line-mesh-based, low-amplitude 
suppression techniques to eliminate spiral-wave turbulence in 
these models for cardiac tissue (Sinha et al, 2001; Pandit et al., 
2002; Shajahan et al, 2009; Majumder et al, 2011a). 

We have first considered simulations in cable-type domains, 
which are ideally suited for the calculation of 0 and X. We find 
that PD leads to a periodic, spatial modulation of 6 and a tempo- 
rally periodic modulation of X; the degrees of these modulations 
depend on the amplitude A x and frequency f x of the PD. To the 



best of our knowledge, such modulations have not been quan- 
tified in any earlier study, although a few studies (Zhang et al., 
2006; Chen et al., 2008) have suggested, in the context of spiral 
waves, that such modulations can arise because of a Doppler-type 
effect (Fenton et al, 2002). 

We have considered three types of initial spiral-wave config- 
urations, IC1, IC2, and 7C3, for both the TP06 and TNNP04 
models. In the TP06 model, these configurations evolve, respec- 
tively, to (a) an RSC state, (b) an SMST state, and (c) an MST 
state, in the absence of PD; in the TNNP04 model they evolve, 
respectively, to (a) an RSC state, (b) an RSZ state, and (c) an 
MST state, in the absence of PD. We have used such initial con- 
ditions because various experimental and computational studies 
(Damle et al, 1992; Davidenko et al, 1992; Pertsov et al, 1993; 
Gray et al, 1995; Gray and Jalife, 1996; Beaumont et al, 1998; 
Gray et al, 1998; Witkowski et al., 1998; Qu et al, 2000a; Fenton 
et al., 2002; Ten Tusscher and Panfilov, 2006; Shajahan et al, 2009) 
have suggested that spiral-wave dynamics in cardiac tissue can 
lead to (a) a stable rotor (Davidenko et al., 1992; Pertsov et al., 
1993; Beaumont et al., 1998), as in our RSC state, (b) a single, 
meandering rotor whose time series is chaotic (Gray et al., 1995; 
Gray and Jalife, 1996; Qu et al., 2000a), as in our SMST state, and 
(c) multiple rotors, which yield a state with spatiotemporal chaos 
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(Damle et al, 1992; Gray et al., 1998; Witkowski et al, 1998; Qu 
et al., 2000a), as in our MST state. Thus, our initial conditions, 
IC1, IC2, and 7C3, lead to the three major types of spiral-wave 
evolutions, and slight variants thereof (e.g., RSZ), which have 
been seen in earlier studies and whose evolution we study now 
with PD. This shows that spiral-wave dynamics in these two mod- 
els, without PD, depends sensitively on the ionic details of these 
models and on the initial conditions. A rich variety of spiral- 
wave behaviors result when we add PD to the TP06 and TNNP04 
models, which are quantitatively different from, but qualitatively 
similar to, that in the TP06 model. These differences arise because 
of the differences in the calcium-ion dynamics in these models; 
such dynamics can play an important role in the spatiotempo- 
ral evolution of spiral waves (Weiss et al., 2005; Ter Keurs and 
Boyden, 2007). Our principal findings here can be summarized 
as follows: In the presence of PD, an RS state may show (a) 
periodic behavior with high-order cycles in time, (b) temporally 
quasiperiodic (QP) evolution, (c) a state with spiral-wave tur- 
bulence, or (d) a quiescent state Q. For an ST state, which can 
be of SMST or MST types, PD can either leave the system in 
an ST state or make it evolve to a Q state, in which all spirals 
either annihilate each other or are absorbed at the boundaries 
of the simulation domain. Precisely which one of these states 
is obtained depends sensitively on our initial conditions and on 
the PD parameters A x , A y , f x , and f y . This sensitive dependence 
on parameters is a hallmark of extended dynamical systems that 
show spatiotemporal chaos (Shajahan et al, 2007, 2009). Thus, 
our study systematizes the effects of PD on spiral-wave dynam- 
ics and turbulence in two, biophysically realistic mathematical 
models for cardiac tissue; and it complements earlier studies of 
spiral-wave dynamics, in such models, that have concentrated on 
the dependence of such dynamics on ion-channel and electro- 
physiological properties (Karma, 1994; Qu et al, 1999, 2000a,b) 
and on conduction (Xie et al., 1998; Ten Tusscher and Panfilov, 
2003; Shajahan et al, 2007, 2009) and ionic inhomogeneities 
(Shajahan et al, 2009; Majumder et al., 2011a,b). By using the 
biophysically realistic TP06 and TNNP04 models for cardiac tis- 
sue, our study generalizes the work of Zhang et al. (2004, 2006) 
and Chen et al. (2008) on spiral-wave instabilities in a simple, 
two- variable model for cardiac tissue, which is subject to PD. 

Moreover, our studies have used three types of spiral-wave 
initial configurations to examine, via extensive and systematic 
numerical calculations, the transitions between different states of 
our system as the amplitude and frequency of the PD are var- 
ied. Our work extends significantly earlier studies of PD (Zhang 
et al, 2004, 2006; Chen et al, 2008) and mechanical deforma- 
tion (Panfilov et al, 2005, 2007; Weise et al., 2011). Zhang et al. 
(2004) have studied the instability of a spiral wave of electrical 
activation by introducing, in a simple, two-variable, FitzHugh- 
Nagumo-type model (Aliev and Panfilov, 1996) for cardiac tissue, 
the possibility of periodic, temporal oscillations in the diffusion 
constant. Their study shows that the resulting deformation can 
lead to a transition from a stable, RS state to an ST state with 
multiple spirals. In another study Zhang et al. (2006) have shown 
that such an ST state can be driven into a quiescent state with 
no spirals, if the oscillation frequency of the PD is chosen to be 
close to the characteristic frequency of the spiral wave in the RS 



state of the system. Chen et al. (2008) have studied the effects 
of PD in the two-variable, Bar model (Bar and Eiswirth, 1993) 
on spiral- wave dynamics by varying the parameter e, which sets 
the time scale of the slow variable in this model, and the ampli- 
tude A and frequency / of the PD; their study shows that the 
RS-ST transition can be effected by changing e, A, and/ suitably. 
Panfilov et al. (2007) have shown that mechanical deformation 
can either (a) induce a spiral wave to drift or (b) break up spi- 
ral waves and thus lead to complex spatiotemporal patterns in the 
three-variable Fenton-Karma model (Fenton and Karma, 1998) 
for cardiac tissue; the model that (Panfilov et al., 2007) use for 
deformation is different from, and more realistic than, the one 
used in Zhang et al. (2004, 2006) and Chen et al. (2008) in so 
far as it includes a stretch-activated current, which accounts for 
the mechano-electrical feedback in cardiac tissue, whose stress 
tensor controls the deformation; their study shows that rotating 
spirals become unstable both because of the stretch-activated cur- 
rent and the deformation of the tissue. In a related study, which 
also includes stretch -activated currents, Weise et al. (2011) have 
shown that such deformation can lead to pacemaker activity in 
a discrete version of the two-variable, Aliev-Panfilov, reaction- 
diffusion model (Aliev and Panfilov, 1996). Note that the studies 
in Zhang et al. (2004, 2006), Panfilov et al. (2007), Chen et al. 
(2008), and Weise et al. (2011) have used only a particular type 
of spiral-wave configuration in their two-variable models, with 
either PD or mechanical deformation. Moreover, Zhang et al. 
(2004, 2006), Panfilov et al. (2005, 2007), Chen et al. (2008), and 
Weise et al. (2011) have focused on simple, two-variable, math- 
ematical models for cardiac tissue. Thus, these studies cannot 
address spiral-wave dynamics in such tissue at the detailed ionic 
level we consider in our work by using state-of-the-art, ionically 
realistic mathematical models for ventricular tissue. Furthermore, 
the authors of Zhang et al. (2004, 2006) and Chen et al. (2008) 
study the effects of PD on spiral-wave dynamics for a limited 
set of initial conditions. For example, in Zhang et al. (2004), the 
authors have studied the behavior of a single rotating spiral (RS) 
in the presence of PD; the authors of Zhang et al. (2006) have 
used a broken-spiral state as an initial configuration to study the 
elimination of spirals from the system in the presence of PD; in 
Chen et al. (2008), the authors have studied an RS initial config- 
uration and its spatiotemporal evolution with PD. The authors of 
Panfilov et al. (2005, 2007) and Weise et al. (2011) have used an 
RS initial configuration to examine the effect of mechanical defor- 
mation on this RS state; they have not investigated the transitions 
between different spiral-wave states. None of these studies have 
carried out the detailed numerical investigations of spiral-wave 
dynamics that we present in our work, which considers a variety 
of initial conditions. 

Furthermore, we have shown that square- and line-mesh- 
based, low-amplitude suppression schemes eliminate spiral-wave 
turbulence in both the TP06 and TNNP04 models in the 
absence of PD; this line-based suppression scheme is a signifi- 
cant improvement over the square-mesh suppression scheme of 
Sinha et al. (2001), Pandit et al. (2002), Shajahan et al. (2009), 
and Majumder et al. (2011a) because it has fewer control lines. 
However, we have found that the line-based scheme works with 
PD only if the PD is applied along one spatial direction. We have 



www.frontiersin.org 



June 2014 | Volume 5 | Article 207 | 13 



Nayak and Pandit 



Spiral-wave dynamics 



then shown that a minor modification of our line-based suppres- 
sion scheme can suppress spiral-wave turbulence: in particular, 
we introduce a rectangular-mesh-based suppression scheme, in 
which we add a few control lines perpendicular to the parallel 
lines of the line-based suppression scheme; this rectangular-mesh 
scheme is also a significant improvement because it uses fewer 
control lines than the one based on a square mesh. We would 
like to emphasize here that no spiral-wave suppression scheme 
has been tried in the presence of PD hitherto; our study is the first 
to address this issue. 

The formation of patterns in reaction-diffusion type system 
with various types of flows have been investigated in Biktashev 
et al. (1998), Leconte et al. (2003), Kuptsov et al. (2005), 
and Ermakova et al. (2009); in particular, break up of spi- 
ral excitation waves has been observed in a moving excitable 
medium as suggested in Biktashev et al. (1998); these studies 
have shown that linear shear flow can cause spiral-wave breakup 
in an excitable medium. Recent studies in Yoshida (2010) and 
Yashin et al. (2012a,b) have investigated pattern formation in 
a gel medium that can be oscillated mechanically; the effect 
of these mechanical oscillations on the underlying spatiotem- 
poral chemical oscillations, because of a Belousov-Zhabotinsky 
(BZ) reaction, can be studied in such systems. Our work here 
presents the cardiac-tissue analogs of such chemical-oscillation 
studies. 

Our study has used methods of non-linear-dynamics to eval- 
uate the effects of periodic deformation on waves of excitation, 
both plane waves and spiral waves, and has obtained there- 
from important physiological implications. In particular, we have 
shown, for the first time, how spiral-wave turbulence can be 
suppressed even in the presence of such deformation. This has 
potential applications in defibrillation because spiral-wave tur- 
bulence in our mathematical models is the analog of ventricular 
fibrillation (Davidenko et al, 1992; Pertsov et al., 1993; Gray et al., 
1998; Jalife et al., 1998). Furthermore, our study has immediate 
implications for experiments on cell cultures with cardiac cells. 

Various stretching devices have been developed to control the 
contraction and expansion of a cell (Barbee et al, 1994; Huang 
et al., 2010; Wang et al., 2010) and a layer of cells in culture (Lee 
et al, 1996; Sotoudeh et al, 1998; Waters et al, 2001; Winter 
et al., 2002; Rana et al., 2008). In these devices, both uniax- 
ial and biaxial (Pang, 2009) stretching methods can be used to 
deform substrates; moreover, in some of these devices, the stretch- 
ing can be applied in a cyclic manner at different frequencies. This 
stretching-induced deformation of the substrate leads, in turn, 
to the deformation of a layer of culture cells that are attached to 
the substrate. Examples of such studies include the following: The 
study in Barbee et al. ( 1 994) has measured the strain that develops 
in cultured vascular smooth muscle cells when they are deformed 
by the stretching of a substrate to which they adhere. The authors 
of Lee et al. (1996) have used a device, which applies homoge- 
neous, equibiaxial strains of 0-10% to a cell-culture substrate, 
to verify quantitatively the transmission of substrate deformation 
to a 2D sheet of cultured cardiac cells. The studies in Sotoudeh 
et al. (1998) have used endothelial cells in tissue culture, on a 
silicon elastic membrane, and have designed an apparatus that 
allows for the control of the magnitude and frequency of the 



dynamical stretching that is applied uniformly to these cells to 
produce equibiaxial dynamical stretches, with area changes rang- 
ing from 0% to 55% and frequencies ranging from 0 to 2 Hz. The 
authors of Waters et al. (2001) have developed a system for the 
imposition of cyclic biaxial strain to stretch cultured pulmonary 
epithelial cells; similar techniques have been used in Winter et al. 
(2002) to study the effects of strain in cell cultures and in vitro 
experiments. The authors of Rana et al. (2008) have studied the 
response to such stretching in cultured neonatal rat atrial car- 
diomyocytes by using a device that can impose homogeneous 
equibiaxial deformation. Other recent studies include those of 
Huang et al. (2010) and Wang et al. (2010), which have stud- 
ied the mechanical activities of living cell, fiber, and tissue by 
applying both equiaxial and uniaxial deformation, and recording 
the dynamics of the response of these systems by using high- 
resolution imaging techniques; the former experiment has used 
fibroblasts and the latter endothelial cells in culture. We suggest 
that such experimental studies of the responses of cell cultures to 
an applied stress can be easily generalized to study the types of 
problems we have concentrated on here. In particular, by impos- 
ing a periodic deformation on cardiac tissue or cell cultures, 
experiments should be able to verify the predictions we have 
made, on the basis of our in silico studies, about the modulations 
of 9 and k in the presence of PD and the effects of PD on spiral- 
wave dynamics, which we have discussed in detail in the previous 
Section. 

We end with some limitations of our model. The first limi- 
tation is that our model does not include stretch-activated cur- 
rents (Panfilov et al, 2005, 2007; Weise et al, 2011) explicitly; 
but V m depends on PD and all ionic currents depend on V m , 
so PD affects all such currents implicitly. Next, the PD in our 
model affects the electrical activation of our medium but it is 
not, in turn, affected by this activation; by contrast, the model 
for mechanical deformation used in Panfilov et al. (2005, 2007) 
and Weise et al. (201 1) allows for electrical feedback to affect such 
deformation; our model does not include soft-tissue mechanics, 
which can be incorporated in mathematical models for cardiac 
tissue by including stress and strain tensors, from elasticity theory 
(Nash and Hunter, 2000; Nash and Panfilov, 2004; Keldermann 
et al, 2009), as in the studies of Panfilov et al. (2005, 2007) 
and Weise et al. (2011); however, these studies use only a two- 
variable model for cardiac tissue and not the ionically realistic 
TP06 or TNNP04 models that we employ. Moreover, because of 
the absence of detailed ion-channel dynamics, the simple, two- 
variable models for cardiac tissue, which have been used in the 
studies of Panfilov et al. (2005, 2007) and Weise et al. (2011), do 
not account for the effects of deformation on ion-channel activ- 
ity and the intracellular calcium concentration as suggested in 
Cherubini et al. (2008), Pathmanathan and Whiteley (2009), and 
Ambrosi et al. (2011). In spite of the simplicity of our model for 
PD, our study captures various features of spiral-wave dynamics 
that have been observed in models that include stretch -activated 
currents (Panfilov et al., 2007); in particular, our model dis- 
plays spiral-wave breakup because of PD. Thus, this result of 
ours is robust. The only qualitative effect that our study misses 
is deformation-induced pacemaker activity, for which it has been 
argued (Kohl et al, 1999; Panfilov et al, 2005; Weise et al, 2011) 
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that stretch-activated currents are essential. To the best of our 
knowledge, our elucidation of the effects of PD on spiral-wave 
dynamics in mathematical models for cardiac tissue, although 
simple in its modeling of PD, is the first study that explores the 
effects of PD on spiral-wave dynamics in ionically realistic math- 
ematical models for ventricular tissue. A complete study of a 
realistic model for deformation, with stretch -activated currents, 
and such ionically realistic mathematical models lies beyond the 
scope of the present paper. While this paper was being prepared 
for publication, a new study appeared on a human-ventricular 
mathematical model with mechanical deformation (Weise and 
Panfilov, 2013). The authors of this study have used the same 
type of deformation that they have employed in their previous 
investigations of mechanical deformation (Weise et al., 2011); 
they have studied spiral-wave dynamics in the context of the 
drifting of a spiral wave as a function of G s (here, G s is the 
maximum conductance of the stretch -activated current). In their 
numerical studies they have found that, in a constantly stretched 
medium, there is an increase of the core size and period of a spi- 
ral wave, but no change in its rotational dynamics; in contrast, 
in the dynamically stretched medium, they observe spiral drift. 
We have found such behaviors, to some extent, in our period- 
ically deformed medium; e.g., the initial condition ICl, which 
leads to a stable rotating spiral with a circular tip trajectory, in 
the absence of PD, can lead to a variety of behaviors, which 
can include slow spiral drift in the presence of PD. Our model 
does not include mechano-electrical feedback in a realistic way, 
as we have described above. Therefore, we have not attempted 
to study how different mechanical stimuli, other than the PD 
we consider, initiate or affect spiral-waves in our model; studies 
of other mechanical stimuli lie beyond the scope of our paper. 
In Weise et al. (2011) it has been noted that both electrical and 
mechanical stimuli can cause the formation of a pacemaker in 
cardiac tissue; and mechanical stimuli can translate the mechan- 
ical energy into an electrical stimulus, as argued in Janse et al. 
(2003) and Cooper et al. (2006). We use a monodomain descrip- 
tion for cardiac tissue; and we do not use an anatomically realistic 
simulation domain (Panfilov and Keener, 1995; Trayanova and 
Tice, 2009), muscle-fiber orientation, and transmural heterogene- 
ity (Majumder et al, 201 lb, 2012); the inclusion of these features 
lies beyond the scope of this study. We note, however, that recent 
studies (Potse et al., 2006) have compared potentials resulting 
from normal depolarization and repolarization in a bidomain 
model with those of a monodomain model; these studies have 
shown that the differences between results obtained from a mon- 
odomain model and those obtained from a bidomain model are 
extremely small. 
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SUPPLEMENTARY MATERIAL 
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Video SI | Spatiotemporal evolution of plane waves in cable-type 
domains, with PD along the axial-direction of the cables, for the TP06 
model, shown via pseudocolor plots of the transmembrane potential V m ; 
the time evolution is shown for 0s < t < 2s; we use 10 frames per second 
(fps); in real time each frame is separated from the succeeding frame by 
8 ms. 

Video S2 | Spiral-wave dynamics in the TP06 and TNNP04 models in the 
absence of PD; the time evolution of pseudocolor plots of the 
transmembrane potential V m is shown for 0 s < t < 4 s; we use 10 frames 
per second (fps); in real time each frame is separated from the succeeding 
frame by 8 ms. 

Video S3 | Spiral-wave dynamics for the TP06 model in the presence of PD 
along both x and y directions with an initial condition of type ICl; the 
time evolution of pseudocolor plots of the transmembrane potential V m is 
shown for 0 s < t < 4 s; we use 10 frames per second (fps); in real time 
each frame is separated from the succeeding frame by 8 ms. 

Video S4 | Spiral-wave dynamics for the TP06 model in the presence of PD 
along both x and y directions with an initial condition of type IC2; the 
time evolution of pseudocolor plots of the transmembrane potential V m is 
shown for 0 s < t < 4 s; we use 10 frames per second (fps); in real time 
each frame is separated from the succeeding frame by 8 ms. 

Video S5 | Spiral-wave dynamics for the TP06 model in the presence of PD 
along both x and y directions with an initial condition of type IC3; the 
time evolution of pseudocolor plots of the transmembrane potential V m is 
shown for 0 s < t < 4 s; we use 10 frames per second (fps); in real time 
each frame is separated from the succeeding frame by 8 ms. 

Video S6 | Spiral-wave dynamics for the TNNP04 model in the presence of 
PD along both x and y directions with an initial condition of type ICl; the 
time evolution of pseudocolor plots of the transmembrane potential V m is 
shown for 0 s < t < 4, s; we use 10 frames per second (fps); in real time 
each frame is separated from the succeeding frame by 8 ms. 

Video S7 | Spiral-wave dynamics for the TNNP04 model in the presence of 
PD along both x and y directions with an initial condition of type IC2; the 
time evolution of pseudocolor plots of the transmembrane potential V m is 
shown for 0 s < t < 4 s; we use 10 frames per second (fps); in real time 
each frame is separated from the succeeding frame by 8 ms. 

Video S8 | Spiral-wave dynamics for the TNNP04 model in the presence of 
PD along both x and y directions with an initial condition of type IC3; the 
time evolution of pseudocolor plots of the transmembrane potential V m is 
shown for 0 s < t < 4 s; we use 10 frames per second (fps); in real time 
each frame is separated from the succeeding frame by 8 ms. 

Video S9 | Comparison of spiral-wave suppression by low-amplitude 
pulses on square, line, and rectangular control meshes in the TP06 model, 
with PD along both x and y directions, for the time interval 0 s < t < 1 s; 
we depict pseudocolor plots of V m and use 10 frames per second (fps); in 
real time each frame is separated from the succeeding frame by 8 ms. 

REFERENCES 

Aliev, R. R., and Panfilov, A. V. (1996). A simple two -variable model of car- 
diac excitation. Chaos Solitons Fractals 7, 293-301. doi: 10.1016/0960-0779(95) 
00089-5 

Ambrosi, D., Arioli, G., Nobile, R, and Quarteroni, A. (2011). Electromechanical 
coupling in cardiac dynamics: the active strain approach. SIAM J. Appl. Math. 
71, 605-621. doi: 10.1137/100788379 

Bar, M., and Eiswirth, M. (1993). Turbulence due to spiral breakup in a con- 
tinuous excitable medium. Phys. Rev. E 48, R1635. doi: 10.1103/PhysRevE.48. 
R1635 

Barbee, K. A., Macarak, E. J., and Thibault, L. E. (1994). Strain measurements in 
cultured vascular smooth muscle cells subjected to mechanical deformation. 
Annal. Biomed. Eng. 22, 14-22. doi: 10.1007/BF02368218 



www.frontiersin.org 



June 2014 | Volume 5 | Article 207 | 15 



Nayak and Pandit 



Spiral-wave dynamics 



Barkley, D., Kness, M., and Tuckerman, L. S. (1990). Spiral-wave dynamics in a sim- 
ple model of excitable media: the transition from simple to compound rotation. 
Phys. Rev. A 42, 2489-2492. doi: 10.1103/PhysRevA.42.2489 

Beaumont, J., Davidenko, N., Davidenko, J. M., and Jalife, J. (1998). Spiral waves in 
two-dimensional models of ventricular muscle: formation of a stationary core. 
Biophys. J. 75, 1-14. doi: 10.1016/S0006-3495(98)77490-9 

Bernus, O., Wilders, R., Zemlin, C. W., Verschelde, H., and Panfilov, A. V. (2002). 
A computationally efficient electrophysiological model of human ventricular 
cells. Am. ]. Physiol. Heart Circ. Physiol. 51, H2296-H2308. doi: 10.1152/ajp- 
heart.0073 1.2001 

Biktashev, V., Holden, A., Tsyganov, M., Brindley, J., and Hill, N. (1998). Excitation 

wave breaking in excitable media with linear shear flow. Phys. Rev. Lett. 81, 2815. 

doi: 10.1103/PhysRevLett.81.2815 
Chen, J.-X., Xu, J.-R., Yuan, X.-R, and Ying, H.-P. (2008). Influences of periodic 

mechanical deformation on spiral breakup in excitable media. /. Phys. Chem. B 

113, 849-853. doi: 10.1021/jp806811j 
Chen, P.-S., Wu, T.-J., Ting, C.-T., Karagueuzian, H. S., Garfinkel, A., Lin, S.- 

E, et al. (2003). A tale of two fibrillations. Circulation 108, 2298-2303. doi: 

10.1161/01 .CIR.0000094404.26004.07 
Chen, S.-A., Hsieh, M.-H., Tai, C.-T., Tsai, C.-E, Prakash, V., Yu, W.-C, et al. 

(1999) . Initiation of atrial fibrillation by ectopic beats originating from the pul- 
monary veins electrophysiological characteristics, pharmacological responses, 
and effects of radiofrequency ablation. Circulation 100, 1879-1886. doi: 
10.1161/01.CIR.100.18.1879 

Cherry, E„ and Fenton, F. (2008). Visualization of spiral and scroll waves in 
simulated and experimental cardiac tissue. New J. Phys. 10, 125016. doi: 
10.1088/1367-2630/10/12/125016 

Cherubini, C, Filippi, S., Nardinocchi, P., and Teresi, L. (2008). An electromechan- 
ical model of cardiac tissue: constitutive issues and electrophysiological effects. 
Prog. Biophys. Mol. Biol. 97, 562-573. doi: 10.1016/j.pbiomolbio.2008.02.001 

Clayton, R., and Panfilov, A. (2008). A guide to modelling cardiac electrical activ- 
ity in anatomically detailed ventricles. Prog. Biophys. Mol. Biol. 96, 19-43. doi: 
10.1016/j.pbiomolbio.2007.07.004 

Cooper, P. J., Epstein, A., MacLeod, I. A., Schaaf, S., Sheldon, J., Boulin, C, et al. 
(2006). Soft tissue impact characterisation kit (stick) for ex situ investigation 
of heart rhythm responses to acute mechanical stimulation. Prog. Biophys. Mol. 
Biol 90, 444-468. doi: 10.1016/j.pbiomolbio.2005.07.004 

Damle, R. S., Kanaan, N. M., Robinson, N. S., Ge, Y.-Z., Goldberger, J. J., 
and Kadish, A. H. (1992). Spatial and temporal linking of epicardial 
activation directions during ventricular fibrillation in dogs. Evidence for 
underlying organization. Circulation 86, 1547-1558. doi: 10. 1161/01. CIR. 
86.5.1547 

Davidenko, J. M., Kent, P. E, Chialvo, D. R., Michaels, D. C, and Jalife, J. (1990). 

Sustained vortex-like waves in normal isolated ventricular muscle. Proc. Natl. 

Acad. Sci. U.S.A. 87, 8785-8789. doi: 10.1073/pnas.87.22.8785 
Davidenko, J. M., Pertsov, A. V., Salomonsz, R., Baxter, W., and Jalife, J. (1992). 

Stationary and drifting spiral waves of excitation in isolated cardiac muscle. 

Nature555, 349-351. doi: 10.1038/355349a0 
De Bakker, J., Van Capelle, E, Janse, M. J., Tasseron, S., Vermeulen, J. T, De Jonge, 

N., et al. (1993). Slow conduction in the infarcted human heart. 'zigzag'course 

of activation. Circulation 88, 915-926. doi: 10.1161/01.CIR.88.3.915 
Ermakova, E. A., Shnol, E. E., Panteleev, M. A., Butylin, A. A., Volpert, V., and 

Ataullakhanov, F. I. (2009). On propagation of excitation waves in moving 

media: the fitzhugh-nagumo model. PLoS ONE 4:e4454. doi: 10.1371/jour- 

nal.pone.0004454 

Fenton, E, and Karma, A. (1998). Vortex dynamics in three-dimensional continu- 
ous myocardium with fiber rotation: filament instability and fibrillation. Chaos 
8, 20-47. doi: 10.1063/1.166311 

Fenton, F. EL, Cherry, E. M., and Glass, L. (2008). Cardiac arrhythmia. Scholarpedia 
3, 1665. doi: 10.4249/scholarpedia.l665 

Fenton, F. H., Cherry, E. M., Hastings, H. M., and Evans, S. J. (2002). Multiple 
mechanisms of spiral wave breakup in a model of cardiac electrical activity. 
Chaos 12, 852-892. doi: 10.1063/1.1504242 

Fogoros, R. N. (2011). Ventricular fibrillation. Heart Dis. Available online at: http:// 
heartdisease. about. com/od/palpitationsarrhythmias/a/Ventricular- Fibrillation, 
htm 

Garfinkel, A, Kim, Y.-H., Voroshilovsky, O., Qu, Z., Kil, J. R., Lee, M.-H., et al. 

(2000) . Preventing ventricular fibrillation by flattening cardiac restitution. Proc. 
Natl. Acad. Sci. U.S.A. 97, 6061-6066. doi: 10.1073/pnas.090492697 



Gray, R., and Jalife, J. (1996). Self-organized drifting spiral waves as a mechanism 

for ventricular fibrillation. Circulation 94, 1-48. 
Gray, R. A., Jalife, J., Panfilov, A. V, Baxter, W. T, Cabo, C, Davidenko, J. M., 

et al. (1995). Mechanisms of cardiac fibrillation. Science 270, 1222-1225. doi: 

10.1126/science.270.5239.1222 
Gray, R. A., Pertsov, A. M., and Jalife, J. (1998). Spatial and temporal organization 

during cardiac fibrillation. Nature 392, 75-78. doi: 10.1038/33537 
Gray, R. A., Wikswo, J. P., and Otani, N. F. (2009). Origin choice and petal loss 

in the flower garden of spiral wave tip trajectories. Chaos 19, 033118. doi: 

10.1063/1.3204256 

Huang, L., Mathieu, P. S., and Helmke, B. P. (2010). A stretching device for 
high-resolution live-cell imaging. Annal. Biomed. Eng. 38, 1728-1740. doi: 
10.1007/sl0439-010-9968-7 

Ideker, R. E., and Rogers, J. M. (2006). Human ventricular fibrillation wan- 
dering wavelets, mother rotors, or both? Circulation 114, 530-532. doi: 
10.1161/CIRCULATIONAHA.106.644765 

Ikeda, T., Uchida, T, Hough, D., Lee, J. J., Fishbein, M. C, Mandel, W. J., et al. 
(1996). Mechanism of spontaneous termination of functional reentry in isolated 
canine right atrium evidence for the presence of an excitable but nonexcited 
core. Circulation 94, 1962-1973. doi: 10.1161/01.CIR.94.8.1962 

Jalife, J., Gray, R. A., Morley, G. E., and Davidenko, J. M. (1998). Self-organization 
and the dynamical nature of ventricular fibrillation. Chaos 8, 79-93. doi: 
10.1063/1.166289 

Janse, M. J., Coronel, R., Wilms-Schopman, F. J., and de Groot, J. R. (2003). 

Mechanical effects on arrhythmogenesis: from pipette to patient. Prog. Biophys. 

Mol Biol 82, 187-195. doi: 10.1016/S0079-6107(03)00015-4 
Karma, A. (1994). Electrical alternans and spiral wave breakup in cardiac tissue. 

Chaos 4, 461-472. doi: 10.1063/1.166024 
Keldermann, R., Nash, M., and Panfilov, A. (2009). Modeling cardiac mechano- 

electrical feedback using reaction-diffusion-mechanics systems. Physica D 238, 

1000-1007. doi: 10.1016/j.physd.2008.08.017 
Keldermann, R. H., Nash, M. P., Gelderblom, H., Wang, V. Y, and Panfilov, A. V. 

(2010). Electromechanical wavebreak in a model of the human left ventricle. 

Am. J. Physiol. Cell Physiol. 299, H134. doi: 10.1152/ajpheart.00862.2009 
Kohl, P., Hunter, P., and Noble, D. (1999). Stretch-induced changes in heart rate 

and rhythm: clinical observations, experiments and mathematical models. Prog. 

Biophys. Mol. Biol. 71, 91-138. doi: 10.1016/S0079-6107(98)00038-8 
Koller, M. L., Riccio, M. L., and Gilmour, R. F Jr. (1998). Dynamic restitution of 

action potential duration during electrical alternans and ventricular fibrillation. 

Am. J. Physiol. Heart Circ. Physiol. 275, H1635-H1642. 
Kuijpers, N. H., Potse, M., van Dam, P. M., ten Eikelder, H. M., Verheule, S., 

Prinzen, F. W., et al. (2011). Mechanoelectrical coupling enhances initiation 

and affects perpetuation of atrial fibrillation during acute atrial dilation. Heart 

Rhythm 8, 429-436. doi: 10.1016/j.hrthm.2010.11.020 
Kuptsov, P. V, Satnoianu, R. A., and Daniels, P. G. (2005). Pattern formation in a 

two-dimensional reaction-diffusion channel with poiseuille flow. Phys. Rev. Ser. 

E72, 036216. doi: 10.1103/PhysRevE.72.036216 
Leconte, M., Martin, J., Rakotomalala, N., and Salin, D. (2003). Pattern of 

reaction diffusion fronts in laminar flows. Phys. Rev. Lett. 90, 128302. doi: 

10.1103/PhysRevLett.90.128302 
Lee, A., Delhaas, T, Waldman, L. K., Mackenna, D. A., Villarreal, F. J., and 

McCulloch, A. D. (1996). An equibiaxial strain system for cultured cells. Am. 

J. Physiol. Cell Physiol. 271, C1400-C1408. 
Lim, Z. Y, Maskara, B., Aguel, E, Emokpae, R., and Tung, L. (2006). Spiral 

wave attachment to millimeter-sized obstacles. Circulation 1 14, 2113-2121. doi: 

10.1161/CIRCULATIONAHA.105.598631 
Majumder, R., Nayak, A. R., and Pandit, R. (2012). Nonequilibrium arrhythmic 

states and transitions in a mathematical model for diffuse fibrosis in human 

cardiac tissue. PLoS ONE 7:e45040. doi: 10.1371/journal.pone.0045040 
Majumder, R., Nayak, A. R., and Pandit, R. (2011a). "An overview of spiral-and 

scroll-wave dynamics in mathematical models for cardiac tissue," in Heart Rate 

and Rhythm, eds O. N. Tripathi, U. Ravens, and M. C. Sanguinetti (Berlin; 

Heidelberg: Springer), 269-282. 
Majumder, R., Nayak, A. R., Pandit, R., Majumder, R., Nayak, A., Pandit, R., 

et al. (2011b). Scroll-wave dynamics in human cardiac tissue: lessons from a 

mathematical model with inhomogeneities and fiber architecture. PLoS ONE 

6:el8052. doi: 10.1371/journal.pone.0018052 
Masse, S., Downar, E., Chauhan, V., Sevaptsidis, E., and Nanthakumar, K. (2007). 

Ventricular fibrillation in myopathic human hearts: mechanistic insights from 



Frontiers in Physiology | Computational Physiology and Medicine June 2014 | Volume 5 | Article 207 | 16 



Nayak and Pandit 



Spiral-wave dynamics 



in vivo global endocardial and epicardial mapping. Am. J. Physiol. Heart Circ. 
Physiol. 292, H2589-H2597. doi: 10.1152/ajpheart.01336.2006 
McCulloch, A. D., Smaill, B. H., and Hunter, P. J. (1987). Left ventricular epicardial 
deformation in isolated arrested dog heart. Am. }. Physiol. Heart Circ. Physiol. 
252, H233-H241. 

Miguel, A., de la Rubia, F. J., and Ivanov, P. C. (2009). Spiral wave annihilation by 
low-frequency planar fronts in a model of excitable media. EPL (Europhys. Lett.) 
86, 18005. doi: 10.1209/0295-5075/86/18005 

Nash, M. P., and Hunter, P. J. (2000). Computational mechanics of the heart. /. 
Elast. Phys. Sci. Solids 61, 113-141. doi: 10.1023/A:1011084330767 

Nash, M. P., Mourad, A., Clayton, R. H., Sutton, P. M., Bradley, C. P., Hayward, M., 
et al. (2006). Evidence for multiple mechanisms in human ventricular fibrilla- 
tion. Circulation 114, 536-542. doi: 10.1161/CIRCULATIONAHA.105.602870 

Nash, M. P., and Panfilov, A. V. (2004). Electromechanical model of excitable tissue 
to study reentrant cardiac arrhythmias. Prog. Biophys. Mol. Biol. 85, 501-522. 
doi: 10.1016/j.pbiomolbio.2004.01.016 

Nayak, A. R. (2013). Spiral-Wave Dynamics in lonically Realistic Mathematical 
Models for Human Ventricular Tissue. Ph.D. thesis, Indian Institute of Science. 

Nayak, A. R., Shajahan, T., Panfilov, A., and Pandit, R. (2013). Spiral-wave dynam- 
ics in a mathematical model of human ventricular tissue with myocytes and 
fibroblasts. PLoS ONE 8:9. doi: 10.1371/journal.pone.0072950 

Noble, D. (2002). Modelling the heart: insights, failures and progress. Bioessays 24, 
1155-1163. doi: 10.1002/bies.l0186 

Otani, N. F. (2002). A primary mechanism for spiral wave meandering. Chaos 12, 
829-842. doi: 10.1063/1.1503921 

Pandit, R., Pande, A., Sinha, S., and Sen, A. (2002). Spiral turbulence and spa- 
tiotemporal chaos: characterization and control in two excitable media. Physica 
A 306, 211-219. doi: 10.1016/S0378-4371(02)00499-5 

Panfilov, A., Keldermann, R., and Nash, M. (2007). Drift and breakup of spiral 
waves in reaction-diffusion-mechanics systems. Proc. Natl. Acad. Sci. U.S.A. 
104, 7922-7926. doi: 10.1073/pnas.0701895104 

Panfilov, A. v., and Keener, J. (1995). Re-entry in an anatomical model of 
the heart. Chaos Solitons Fractals 5, 681-689. doi: 10.1016/0960-0779(93) 
E0050-L 

Panfilov, A. V., Keldermann, R., and Nash, M. (2005). Self-organized pacemakers 
in a coupled reaction-diffusion-mechanics system. Phys. Rev. Lett. 95, 258104. 
doi: 10.1103/PhysRevLett.95.258104 

Pang, Q. (2009). Design and Development of a Biostretch Apparatus for Tissue 
Engineering. Ph.D. thesis, University of Toronto. 

Pathmanathan, P., and Whiteley, J. P. (2009). A numerical method for car- 
diac mechanoelectric simulations. Annul. Biomed. Eng. 37, 860-873. doi: 
10. 1007/sl0439-009-9663-8 

Pertsov, A. M., Davidenko, J. M., Salomonsz, R., Baxter, W. T., and Jalife, J. (1993). 
Spiral waves of excitation underlie reentrant activity in isolated cardiac muscle. 
Circ. Res. 72, 631-650. doi: 10.1161/01.RES.72.3.631 

Potse, M., Dube, B., Richer, J., Vinet, A., and Gulrajani, R. M. (2006). A comparison 
of monodomain and bidomain reaction-diffusion models for action potential 
propagation in the human heart. IEEE Trans. Biomed. Eng. 53, 2425-2435. doi: 
10. 1 109/TBME.2006.880875 

Qu, Z., Weiss, J. N., and Garfinkel, A. (1999). Cardiac electrical restitution prop- 
erties and stability of reentrant spiral waves: a simulation study. Am. J. Physi. 
Heart Circ. Physiol. 276, H269-H283. 

Qu, Z., Weiss, J. N., and Garfinkel, A. (2000a). From local to global spa- 
tiotemporal chaos in a cardiac tissue model. Phys. Rev. E 61, 727. doi: 
10.1103/PhysRevE.61.727 

Qu, Z., Xie, R, Garfinkel, A., and Weiss, J. N. (2000b). Origins of spiral wave mean- 
der and breakup in a two-dimensional cardiac tissue model. Annul. Biomed. Eng. 
28, 755-771. doi: 10.1114/1.1289474 

Rana, O. R., Zobel, C, Saygili, E., Brixius, K., Gramley, E, Schimpf, T., et al. 
(2008). A simple device to apply equibiaxial strain to cells cultured on flex- 
ible membranes. Am. ]. Physiol. Heart Circ. Physiol. 294, H532-H540. doi: 
10.1152/ajpheart.00649.2007 

Roger, V. L., Go, A. S., Lloyd-Jones, D. M., Adams, R. J., Berry, J. D., Brown, 
T. M., et al. (2011). Heart disease and stroke statistics-2011 update a 
report from the american heart association. Circulution 123, el8-e209. doi: 
10.1161/CIR.0b013e3182009701 

Roger, V. L., Go, A. S., Lloyd-Jones, D. M., Benjamin, E. J., Berry, J. D., 
Borden, W. B., et al. (2012). Heart disease and stroke statistics-2012 update 



a report from the american heart association. Circulution 125, e2-e220. doi: 
10.1161/CIR.0b013e31823ac046 
Samie, F. H., and Jalife, J. (2001). Mechanisms underlying ventricular tachy- 
cardia and its transition to ventricular fibrillation in the structurally 
normal heart. Cardiovasc. Res. 50, 242-250. doi: 10.1016/S0008-6363(00) 
00289-3 

Shajahan, T., Nayak, A. R., and Pandit, R. (2009). Spiral-wave turbulence and its 
control in the presence of inhomogeneities in four mathematical models of 
cardiac tissue. PLoS ONE 4:e4738. doi: 10.1371/journal.pone.0004738 

Shajahan, T., Sinha, S., and Pandit, R. (2007). Spiral-wave dynamics depend sen- 
sitively on inhomogeneities in mathematical models of ventricular tissue. Phys. 
Rev. £ 75, 0 1 1929. doi: 1 0. 1 103/PhysRevE.75.0 1 1 929 

Sinha, S., Pande, A., and Pandit, R. (2001). Defibrillation via the elimination of 
spiral turbulence in a model for ventricular fibrillation. Phys. Rev. Lett. 86, 
3678-3681. doi: 10.1103/PhysRevLett.86.3678 

Sinha, S., and Sridhar, S. (2007). "Controlling spatiotemporal chaos and spiral tur- 
bulence in excitable media," in Handbook of Chaos Control, 2nd Edn., eds E. 
Scholl and H. G. Schuster (Weinheim: Wiley Online Library), 703-718. doi: 
10.1002/9783527622313.ch32 

Sotoudeh, M., Jalali, S., Usami, S., Shyy, J. Y., and Chien, S. (1998). A strain device 
imposing dynamic and uniform equi-biaxial strain to cultured cells. Annal. 
Biomed. Eng. 26, 181-189. doi: 10.1114/1.48 

Sridhar, S., and Sinha, S. (2008). Controlling spatiotemporal chaos in excitable 
media using an array of control points. EPL (Europhys. Lett.) 81, 50002. doi: 
10.1209/0295-5075/81/50002 

Tang, G., Deng, M., Hu, B., and Hu, G. (2008). Active and passive con- 
trol of spiral turbulence in excitable media. Phys. Rev. E 77, 046217. doi: 
10.1 103/PhysRevE.77.0462 1 7 

Ten Tusscher, K., Noble, D., Noble, P., and Panfilov, A. (2004). A model for human 
ventricular tissue. Am. J. Physiol. Heurt Circ. Physiol. 286, H1573-H1589. doi: 
10.11 52/ajpheart.00794.2003 

Ten Tusscher, K., and Panfilov, A. (2003). Influence of nonexcitable cells on spiral 
breakup in two-dimensional and three-dimensional excitable media. Phys. Rev. 
E 68, 062902. doi: 10. 1 103/PhysRevE.68.062902 

Ten Tusscher, K., and Panfilov, A. (2006). Alternans and spiral breakup in a 
human ventricular. Am. J. Physiol. Heurt Circ. Physiol. 291, H1088-H1 100. doi: 
10.1152/ajpheart.00109.2006 

Ter Keurs, H. E., and Boyden, P. A. (2007). Calcium and arrhythmogenesis. Physiol. 
Rev. 87, 457-506. doi: 10.1152/physrev.00011.2006 

Thompson, S. A., Copeland, C. R., Reich, D. H., and Tung, L. (2011). 
Mechanical coupling between myofibroblasts and cardiomyocytes slows elec- 
tric conduction in fibrotic cell monolayers. Circulution 123, 2083-2093. doi: 
10.1161/CIRCULATIONAHA.110.015057 

Trayanova, N. A., and Tice, B. M. (2009). Integrative computational models of car- 
diac arrhythmias-simulating the structurally realistic heart. Drug Discov. Today 
Dis. Models 6, 85-91. doi: 10.1016/j.ddmod.2009.08.001 

Valderrabano, M., Kim, Y.-H, Yashima, M., Wu, T.-J., Karagueuzian, H. S., and 
Chen, P.-S. (2000). Obstacle-induced transition from ventricular fibrillation to 
tachycardia in isolated swine right ventriclesinsights into the transition dynam- 
ics and implications for the critical mass. /. Am. Coll. Cardiol. 36, 2000-2008. 
doi: 10.1016/S0735-1097(00)00941-4 

Wang, D., Xie, Y, Yuan, B., Xu, J., Gong, P., and Jiang, X. (2010). A stretching device 
for imaging real-time molecular dynamics of live cells adhering to elastic mem- 
branes on inverted microscopes during the entire process of the stretch. Integr. 
Biol. 2, 288-293. doi: 10.1039/b920644b 

Waters, C. M., Glucksberg, M. R., Lautenschlager, E. P., Lee, C.-W., Van Matre, 
R. M., Warp, R. J., et al. (2001). A system to impose prescribed homogenous 
strains on cultured cells. /. Appl. Physiol. 91, 1600-1610. 

Weise, L. D., Nash, M. P., and Panfilov, A. V. (2011). A discrete model to study 
reaction-diffusion-mechanics systems. PLoS ONE 6:e21934. doi: 10.1371/jour- 
nal.pone.0021934 

Weise, L. D., and Panfilov, A. V. (2013). A discrete electromechanical model for 
human cardiac tissue: effects of stretch-activated currents and stretch conditions 
on restitution properties and spiral wave dynamics. PLoS ONE 8:e59317. doi: 
10.1371/annotation/9ceadf50-eb8f-4051-9e41-772884d47385 

Weiss, J. N., Qu, Z., Chen, P.-S., Lin, S.-F., Karagueuzian, H. S., Hayashi, H., et al. 
(2005). The dynamics of cardiac fibrillation. Circulation 112, 1232-1240. doi: 
10.1161/CIRCULATIONAHA.104.529545 



www.frontiersin.org 



June 2014 | Volume 5 | Article 207 | 17 



Nayak and Pandit 



Spiral-wave dynamics 



Winter, L. C, Gilbert, J. A., Elder, S. H., and Bumgardner, J. D. (2002). A device for 
imposing cyclic strain to cells growing on implant alloys. Annal. Biotned. Eng. 
30, 1242-1250. doi: 10.1114/1.1529195 

Witkowski, F. X., Leon, L. J., Penkoske, P. A., Giles, W. R., Spano, M. L., Ditto, 
W. L., et al. (1998). Spatiotemporal evolution of ventricular fibrillation. Nature 
392, 78-82. doi: 10.1038/32170 

Wu, T.-J., Lin, S.-E, Baher, A., Qu, Z., Garfmkel, A., Weiss, J. N., et al. (2004). 
Mother rotors and the mechanisms of d600-induced type 2 ventricular fibril- 
lation. Circulation 110, 2110-2118. doi: 10.1161/01.CIR.0000143834.51102.91 

Xie, E, Qu, Z., and Garfinkel, A. (1998). Dynamics of reentry around a circular 
obstacle in cardiac tissue. Phys. Rev. E 58, 6355. doi: 10.1103/PhysRevE.58.6355 

Xie, E, Qu, Z., Weiss, J. N., and Garfinkel, A. (2001). Coexistence of multiple spiral 
waves with independent frequencies in a heterogeneous excitable medium. Phys. 
Rev. £63, 031905. doi: 10.1103/PhysRevE.63.031905 

Yashin, V. V., Kuksenok, O., Dayal, P., and Balazs, A. C. (2012a). Mechano- 
chemical oscillations and waves in reactive gels. Rep. Prog. Phys. 75, 066601. doi: 
10.1088/0034-4885/75/6/066601 

Yashin, V. V., Suzuki, S., Yoshida, R., and Balazs, A. C. (2012b). Controlling the 
dynamic behavior of heterogeneous self-oscillating gels. /. Mater. Chem. 22, 
13625-13636. doi: 10.1039/c2jm32065g 

Yoshida, R. (2010). Self-oscillating gels driven by the belousov-zhabotinsky 
reaction as novel smart materials. Adv. Mater. 22, 3463-3483. doi: 
10.1002/adma.200904075 

Yuan, G, Wang, G, and Chen, S. (2005). Control of spiral waves and spatiotempo- 
ral chaos by periodic perturbation near the boundary. EPL (Europhys. Lett.) 72, 
908. doi: 10.1209/epl/i2004-10553-2 

Zhang, H., Hu, B., andHu, G. (2003). Suppression of spiral waves and spatiotempo- 
ral chaos by generating target waves in excitable media. Phys. Rev. E 68, 026134. 
doi: 10.1103/PhysRevE.68.026134 



Zhang, H., Li, B.-W, Sheng, Z.-M., Cao, Z., and Hu, G. (2006). The effect of 

mechanical deformation on spiral turbulence. EPL (Europhys. Lett.) 76, 1109. 

doi: 10.1209/epl/i2006-10391-2 
Zhang, EL, Ruan, X.-S., Hu, B., and Ouyang, Q. (2004). Spiral breakup due 

to mechanical deformation in excitable media. Phys. Rev. E 70, 016212. doi: 

10.1103/PhysRevE.70.016212 
Zimmermann, M., and Kalusche, D. (2001). Fluctuation in autonomic tone is a 

major determinant of sustained atrial arrhythmias in patients with focal ectopy 

originating from the pulmonary veins. /. Cardiovasc. Electrophysiol. 12, 285-291. 

doi: 10.1046/j.l540-8167.2001.00285.x 
Zipes, D. P., and Wellens, H. J. (1998). Sudden cardiac death. Circulation 98, 

2334-2351. doi: 10.1161/01.CIR.98.21.2334 

Conflict of Interest Statement: The authors declare that the research was con- 
ducted in the absence of any commercial or financial relationships that could be 
construed as a potential conflict of interest. 

Received: 29 September 2013; accepted: 14 May 2014; published online: 10 June 2014. 
Citation: Nayak AR and Pandit R (2014) Spiral-wave dynamics in ionically realistic 
mathematical models for human ventricular tissue: the effects of periodic deformation. 
Front. Physiol. 5:207. doi: 10.3389/fphys.2014.00207 

This article was submitted to Computational Physiology and Medicine, a section of the 
journal Frontiers in Physiology. 

Copyright © 2014 Nayak and Pandit. This is an open-access article distributed under 
the terms of the Creative Commons Attribution License (CC BY). The use, distribu- 
tion or reproduction in other forums is permitted, provided the original author(s) 
or licensor are credited and that the original publication in this journal is cited, in 
accordance with accepted academic practice. No use, distribution or reproduction is 
permitted which does not comply with these terms. 



Frontiers in Physiology | Computational Physiology and Medicine 



June 2014 | Volume 5 | Article 207 | 18 



